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Abstract. The structure of the excited 2'^S and 2'^P triplet states of ^He and *He in an applied magnetic 
field B is studied using different approximations of the atomic Hamiltonian. All optical transitions (line 
positions and intensities) of the 1083 nm 2''S-2^P transition are computed as a function of B. The effect 
of metastability exchange collisions between atoms in the ground state and in the 2'^S metastable state is 
studied, and rate equations are derived, for the populations these states in the general case of an isotopic 
mixture in an arbitrary field B. It is shown that the usual spin-temperature description remains valid. A 
simple optical pumping model based on these rate equations is used to study the B-dependence of the 
population couplings which result from the exchange collisions. Simple spectroscopy measurements are 
performed using a single-frequency laser diode on the 1083 nm transition. The accuracy of frequency scans 
and of measurements of transition intensities is studied. Systematic experimental verifications are made 
for B=0 to 1.5 T. Optical pumping effects resulting from hyperfine decoupling in high field are observed 
to be in good agreement with the predictions of the simple model. Based on adequately chosen absorption 
measurements at 1083 nm, a general optical method to measure the nuclear polarisation of the atoms in 
the ground state in an arbitrary field is described. It is demonstrated at B ~0.1 T, a field for which the 
usual optical methods could not operate. 

PACS. 32.60.-|-i Zeeman and Stark effects - 32.70.-n Intensities and shapes of atomic spectral lines - 
32.80.Bx Level crossing and optical pumping 

1 Introduction be at atmospheric pressure to be inhaled in order to per- 

form lung MRI. Adding a neutral buffer gas after comple- 
Highly polarised ^He is used for several apphcations in of OP is a simple method to increase pressure, but 

various domains, for instance to prepare polarised targets ^'^sults in a large dilution of the polarised hehum. Polari- 
for nuclear physics experiments §, to obtain spin filters nation preserving compression of the helium gas after OP 
for cold neutrons iJl, or to perform magnetic resonance "smg different compressing devices is now performed by 
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in the structures of the different excited levels of helium. 
In order to populate the metastable state and per- 
form OP, a plasma discharge is sustained in the helium 
gas. In the various atomic and molecular excited states 
which are populated in the plasma, hyperfine interaction 
transfers nuclear orientation to electronic spin and orbital 
orientations. The electronic angular momentum is in turn 
converted into circular polarisation of the emitted light 
or otherwise lost during collisions. This process is actu- 
ally put to use in the standard optical detection technique 
P,p^ in which the circular polarisation of a chosen helium 
spectral line emitted by the plasma is measured and the 
nuclear polarisation is inferred. The decoupling effect of 
an applied magnetic field unfortunately reduces the sen- 
sitivity of the standard optical detection method above 
10 mT, and hence a different measurement technique must 
be used in high fields. This transfer of orientation also has 
an adverse effect in OP situations by inducing a net loss 
of nuclear polarisation in the gas. The decoupling effect 
of an applied field reduces this polarisation loss and may 
thus significantly improve the OP performance in situ- 
ations of limited efficiency, such as low temperatures or 
high pressures. At low temperature (4.2 K), a reduced 
metastability exchange cross section sets a tight bottle- 
neck and strongly limits the efficiency of OP [[l8yi9| , |20|] ; a 
field increase from 1 to 40 mT was observed to provide an 
increase in nuclear polarisation from 17 % to 29 % in this 
situation At high gas pressures (above a few mbar) 
the proportion of metastable atoms is reduced and the 
creation of metastable molecular species is enhanced, two 
factors which tend to reduce the efficiency of OP ; it is 
not surprising that a significant improvement is obtained 
by suppressing relaxation channels in high field [l^ . 

A systematic investigation of various processes rele- 
vant for OP in non-standard conditions (high field and/or 
high pressure) has been made, and results will be reported 
elsewhere p^ , p3t . As mentioned earlier, an optical mea- 
surement method of nuclear polarisation in arbitrary field 
must be developed. It is based on absorption measure- 
ments of a probe beam, and requires a detailed knowledge 
of magnetic field effects on the 1083 nm transition lines. 
In this article we report on a detailed study of the effect of 
an applied magnetic field on the structure of the 2^S and 
2^P atomic levels, both in '^He and ^He. In the theoreti- 
cal section we first present results obtained using a simple 
effective Hamiltonian and discuss the accuracy of com- 
puted line positions and intensities at various magnetic 
fields, then discuss the effect of an applied field on the 
metastability exchange collisions between the 2'^S state 
and the l^So ground state level of helium. In the experi- 
mental section we present results of measurements of the 
line positions and intensities up to 1.5 T, then describe an 
optical measurement technique of the nuclear polarisation 
of 3 He. 

Let us finally mention that these studies have been 
principally motivated by OP developments, but that their 
results could be of interest to design or interpret experi- 
ments performed on helium atoms in the 2'^S or 2'^P state 
in an applied field, as long as metrological accuracy is not 



required. This includes laser cooling and Zeeman slowing 
[ p4| of a metastable atom beam, magnetic trapping and 
evaporative cooling which recently allowed two groups to 
obtain Bose Einstein condensation with metastable *He 
atoms p^ , p6| , and similar experiments which may probe 
the infiuence of Fermi statistics with ultracold metastable 
■^He atoms. 



2 Theoretical 

The low-lying energy states of helium are represented in 
figure |l|. The fine-structure splittings of the 2'^P state of 
*He and the additional effect of hyperfine interaction on 
splittings of the 2'^S and 2'^P states of "^He are indicated 
in null magnetic field. The ^He 2'^P fine-structure inter- 
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Fig. 1. Ground state and first excited states of He (not to 
scale) , showing the effects of fine and hyperfine interactions on 
excited level structures. The Zeeman effect is assumed negligi- 
ble, except in the ground state of ^He in which it is grossly ex- 
aggerated to highlight the existence of two magnetic sublevels 
and the possibility of nuclear polarisation. Accurate values of 
all splittings are given in the tables in the appendix. 



vals have been computed and measured with a steadily 
improved accuracy p^p^ , in order to provide a test of 
QED calculations and to measure the fine-structure con- 
stant a. High precision isotope shift measurements of the 
2'^S — 2^P transition of helium have also been performed 
to further test QED calculations and to probe the ^He 
nuclear charge radius with atomic physics experiments 
[ p9| . High precision measurements and calculations of the 
Zeeman effect in the 2'^P state of helium have also been 
performed |3C|,|3l|,|32| . They were required in particular to 
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analyse the results of fine-structure intervals measured us- 
ing microwave transitions in an applied magnetic field in 
■^He [H and ^He ||. 

All this work has produced a wealth of data which may 
now be used to compute the level structures and transi- 
tion probabilities with a very high accuracy up to high 
magnetic fields (several Tesla), in spite of a small per- 
sisting disagreement between theory and experiments for 
the g-factor of the Zeeman effect in the 2"^? state ||3^ , |35| . 
However, such accurate computation can be very time- 
consuming, especially for '^He in which fine-structure, hy- 
perfine and Zeeman interaction terms of comparable im- 
portance have to be considered. A simplified approach was 
proposed by the Yale group |36), based on the use of an 
effective Hamiltonian. The hyperfine mixing of the 
and singlet 2^P states is the only one considered in this 
effective Hamiltonian, and its parameters have been ex- 
perimentally determined . A theoretical calculation of 
the "^He 2'^P structure with an accuracy of order 1 MHz 
later confirmed the validity of this phenomenological ap- 
proach In the present article, we propose a further 
simplification of the effective Hamiltonian introduced by 
the Yale group, in which couplings to the singlet states are 
not explicitly considered. Instead we implicitly take into 
account these coupling terms, at least in part, by using the 
2^P splittings measured in zero field to set the eigenvalues 
of the fine-structure matrices. 

Results obtained using this si mplifie d effective Hamil- 
tonian will be compared in section 2.1.3 to those of the one 
including the couplings to the 2^P levels, with 3 or 6 addi- 
tional magnetic sublevels, depending on the isotope. The 
effect of additional terms in a more elaborate form of the 
Zeeman Hamiltonian than the linear approximation 
which we use will also be evaluated. 



2.1 A simple effective Hamiltonian 

2.1.1 Notations 

The l^S ground level of helium is a singlet spin state 
(S'=0), with no orbital angular momentum (L=0), and 
hence has no total electronic angular momentum (J=0). 
In ^He, the nucleus thus carries the only angular momen- 
tum /, giving rise to the two magnetic sublevels m/=±l/2. 
Their relative populations, (1 ± M) /2, define the nuclear 
polarisation M of the ground state. 

The level structure of '^He excited states is determined 
from the fine-structure term and the Zeeman term in the 
Hamiltonian. The fine-structure term Hf^ is easily ex- 
pressed in the total angular momentum representation (J) 
using parameters given in tables ^ and ^ in the Appendix. 
For the Zeeman term in the applied magnetic field B we 
shall use the simple linear form : 



which the 3 sublevels \ms) (with m5=±l,0) are eigen- 
states at any field. For the 2^P state the Zeeman term 
couples sublevels of different J (which however have the 
same mj) and is more easily expressed in the decoupled 
{L,S) representation. We note |mL,ms) the vectors of 
the decoupled basis used in the {L, S) representation and 
I J; mj) those of the coupled basis, which are the zero-field 
9 eigenstates of the 2'^P state. We call Zi to Zg (by increas- 
ing value of the energy, see figure the 9 sublevels of the 
2'^P state in arbitrary field. To simplify notations used in 
the following, we similarly call Yi to Y3 the 3 sublevels of 
the 2^8 state. 
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Fig. 2. Left : magnetic sublevels of *He in the (bottom) 
and 2^P (top) levels in very low magnetic field. J is thus a good 
quantum number and m,j is used to differentiate the atomic 
states. When a small field is applied, the sublevel degeneracies 
are removed and energies increase for increasing values of mj 
(i.e. from Zi to Zg). Right : magnetic sublevels of ^He in very 
low magnetic field. F and mp are now good quantum numbers 
to differentiate the atomic states. Only the -F=3/2 Fx and Py 
states involve a strong mixing of states of different J values . 
Energies increase for increasing values of mp for all sublevels 
except in the 2'^Po state (i.e. from Ai to Ae and Bi to Bie). 



(1) 



The values of Bohr magneton /i^, and of the orbital and 
spin (^-factors g'j^ and g'g are given in table @ in the Ap- 
pendix. The situation is quite simple for the 2^S state, for 



The transition probabilities of all components of the 
1083 nm line are evaluated using the properties of the elec- 
tric dipole transition operator. For a monochromatic laser 
light with frequency w/27r, polarisation vector e\, and in- 
tensity lias (expressed in W/m^), the photon absorption 
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rate for a transition from to Zj is given by : 



1 47ra/ 



{r/2f 



2 ^ij 



(2) 



where a is the fine-structure constant, / the oscillator 
strength of the 2^8 -2^? transition (/=0.5391 [|8[|3|]), 
the electron mass, r'/2 the total damping rate of the op- 
tical coherence between the 2^S and the 2'^P states, and 
Tj^j\ex) the transition matrix element between Yi and 
Zj 0] for the light polarisation vector e\. The Lorentzian 
factor in equation ^ is responsible for resonant absorption 
(the absorption rate is decreased when the light frequency 



differs from the transition frequency ulf^ j^n, fy^l'^' being 
the energy difference between and Zj). At low enough 
pressure, F' is equal to the radiative decay rate F of the 
2^P state (the inverse of its lifetime) , related to the oscil- 
lator strength / by : 



,(4) 



F = 



(3) 



from which 7^=1.022x10 s^ is obtained. Atomic colli- 
sions contribute to F' by a pressure-dependent amount of 
order 10® s~"'^/mbar ^^l- The transition matrix elements 
T^^\e\) are evaluated in the (L, S) representation for each 
of the three light polarisation states X = x±iy (circularly 
polarised light a± propagating along the 2;-axis set by the 
field B) and X — z (transverse light with tt polarisation) 
using the selection rules given in table in the Appendix. 
The transformation operator P4 given in table ^ in the 
Appendix is used to change between the {L, S) and the 
(J) representations in the 2'^P state of ''He and thus write 
a simple 9x9 matrix expression for the Hamiltonian in the 
{L, S) representation, Hg : 



H9 = P^^H9fsP4 + H9Z, 



(4) 



where ifg/s is the diagonal matrix of table ^ in the Ap- 
pendix and Hqz the diagonal matrix written from equa- 
tion 0. 

The level structure of ^He excited states is determined 
from the total Hamiltonian, including hyperfine interac- 
tions. The fine-structure term Hf^ is expressed in the to- 
tal angular momentum representation (J) using slightly 
different parameters (table ^ in the Appendix) due to the 
small mass dependence of the fine-structure intervals. For 
the Zeeman term we also take into account the nuclear 
contribution: 



H^^t^B ig'LL + g'sS + giI)-B. 



(5) 



The values of all ^-factors for '^He are listed in table H 
in the Appendix. For the hyperfine term, we consider the 
simple contact interaction between the nuclear and elec- 
tronic spins and a correction term H^Ys • 



hfs 



A.I-S 



Tjcor 
^hfs- 



(6) 



The matrix elements of the main contact interaction term 
in the decoupled {S, I) representation and the values of the 
constants As and Ap for the 2'^S and 2'^P states are given 
in table H in the Appendix. The correction term only exists 
for the 2P state |3|,|^. When restricted to the triplet 
levels, it only depends on 2 parameters which are given 
with the matrix elements in table [l^ in the Appendix. 

As in the case of ^He, we compute level structure in 
the decoupled representation, now {L, S, I) . The vectors 
of the decoupled bases are noted \ms : m/) for the 2'^S 
state, and \mL,ms : mj) for the 2'^P state. For the 2^S 
state, the 6x6 matrix expression iJg of the Hamiltonian 



He — FIqz + Fq, (7) 

where Hqz is the diagonal matrix written from equation |^ 
and ^6 the matrix representation of Hh fs (table ^ in the 
Appendix). For the 2'^P state, the 18x18 matrix His is 
constructed from iJg/s, Fq, 



H^^Ys and 



m'j \His\ mL,ms : mj) 



S {m'j,mi) {m'i^,m's \Hgfs\mL,ms) 
+ S {m'^,mL) {m'g : m'j \Fe\ms : mj) 



+ ( rrij^, rug : rrij 



Tjcor 
^hfs 



H 



(3) 



The mass-dependent constants in iJg/s and Fq are set to 
the values given for '^He in the Appendix. 

Notations similar to those of the zero-field calculation 
of reference are used. The 6 sublevels of the 2'^S state 
are called Ai to Ag, the 18 levels of the 2'^P state Bi 
to B18 by increasing values of the energy (see figure H)|^. 
An alternative notation system, given in table |l2 in the 
Appendix, is also used for the 6 sublevels of the 2^S state 
when an explicit reference to the total angular momentum 
projection mp is desired. Transition probabilities from A^ 
to Bj due to monochromatic light are given by a formula 
similar to equation : 



{F'/2y 



1 Anaf ^ 



T^Jiex), (8) 



in which the level energy differences fkoij and the tran- 
sition matrix elements Tij{e\) are generalisations of the 
-6=0 values computed in reference [[7|. To simplify nota- 
tions, no upper index (3) is attached to r^j, cuij and Tij in 
the case of ^He atoms. The small mass effect on the oscilla- 
tor strength / (of order 10"* [^) is neglected. Sum rules 
on the transition matrix elements for each polarisation e\ 
give the relations: 



E 



E 



. T,j (ex) = 6. 



(9) 



^ There is a difference with notations used in reference Q : 
within each level of given F, labelling of state names was in- 
creased for convenience from largest to lower values of rriF. In 
an applied field, this would correspond to an energy decrease 
inside each of the F levels (except the Po, F=l/2 level). This 
is the reason for the new labelling convention, which is more 
convenient in an applied field. 
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In the following, the dependence of the transition ma- 
trix elements on the polarisation vector e\ will not be 
explicitly written. Given the selection rules described in 
the Appendix, the polarisation vector corresponding to a 
non-zero transition matrix element can be unambiguously 
determined. 

All optical transition energies and hjJij will be 

referenced to hio^'ji^), the energy of the Ci transition in 
null magnetic field, which connects levels A5 and Ag to 
levels By to Biq. Energy differences will be noted e, for 
instance : 



'He 



,(4) /fc _,,(4) 



-^57(0), eij/h = - (10) 



They are determined from computed energy splittings in 
the and states of each isotope. The isotope shift 
contribution is adjusted to obtain the precisely measured 
C9-D2 interval (810.599 MHz, pi). 



2.1.2 Numerical Results 

Numerical computation of level structure and absorption 
spectra in an applied field B is performed in 3 steps. Matri- 
ces i?9 (for ■*He), and i/ig (for -^He) are computed first. 
A standard matrix diagonalisation package (the double- 
precision versions of the jacobi and eigsrt routines ^l|) is 
then used to compute and sort all eigenvalues (energies) 
and eigenvectors (components of the atomic states in the 

decoupled bases) . All transition matrix elements T^^"^ and 
Tij are finally evaluated, and results are output to files for 
further use or graphic display. Actual computation time 
is insignificant (e.g. 20 ms per value of i? on a PC using 
a compiled Fortran program |Q), so that we chose not 
to use the matrix symmetries to reduce the matrix sizes 
and the computational load, in contrast with references 
|^,|^. We thus directly obtain all transitions of the 1083 
nm line: 19 transitions for ^He (6 for tT_|_ and for (7_, 7 
for TT light polarisation), and 70 transitions for '^He (22 for 
(T_|_ and for (t_, 26 for tt light polarisation). These num- 
bers are reduced for B=Q to 18 for ^He (the Y2 — > Z7, 
i.e. |0) — + |1; 0) TT transition has a null probability) and to 
64 for ^Hc (the F=l/2^ F=5/2 transitions, from A5 or 
Ag to any of Bi to Bg are forbidden by selection rules). 
Due to level degeneracy, the usual Dq, Diand D2 lines of 
''He, and Ci to Cg lines of '^He are obtained for B=Q 
(figure ||, and table ^ in the Appendix) . 

An applied magnetic field removes level degeneracies 
(Zeeman splittings appear) and modifies the atomic states 
and hence the optical transition probabilities. Examples 
of level Zeeman energy shifts are shown in figure for 
the 2^S state and the highest-lying 2'^P levels (originating 
from the 2'^Po state at low field). In the 2^S state, the 
two Zeeman splittings are simply proportional to B (28 
GHz/T) for ^He, as inferred from equation ll|. For ■^He, 
significant deviations from linearity occur for A3 and 
A5 above a field of order 0.1 T, for which Zeeman shifts 
are significant compared to the hyperfine splitting. Level 
crossing of A4 and A5 (and hence interchange of names 
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and ^He (right) for B=0 



GHz 

Components of the 1083 nm line in ''He (left) 
Dj is the usual name of the tran- 
sition to the 2'^P,j level. Ci to Cg refer to the transitions by 
increasing energy. Co is the additional lowest energy compo- 
nent corresponding to a forbidden transition in zero field (see 
text). Bottom : positions of all the spectral lines resulting from 
level splittings and isotope shift. The total frequency span from 
Ci to Do is 70.442 GHz. 



of the two states) occurs at 0.1619 T. A high-field decou- 
pled regime is almost reached for 1.5 T (figure^, bottom). 
Analytical expression for state energies and components 
on the decoupled {\ms : toj)} basis are given in the Ap- 
pendix (equations |A10 to A17). In the 2'^P state, the sit- 
uation is more complex due to fine-structure interactions 
and to the larger number of levels. In particular no level 
remains unaffected by B in contrast with the situation for 
Y2. As a result, tt transitions in ''He from Y2 to Z3 or 
Zg have non-zero Zeeman frequency shifts, which are pro- 
portional to B^ at low and moderate field although they 
originate from the linear Zeeman term in the Hamiltonian 
(equation |l|). Even the 2'^Po states, for which level mixing 
is weakest due to the large fine-structure gap, experience 
significant field effects (figure ^, top and top insert). The 
spHtting between B17 and Big (2.468 GHz/T at low field) 
is only weakly affected by the nuclear term in equation H 
(±16.2 MHz/T), and mostly results from level mixing and 
electronic Zeeman effects. 



Examples of optical transition intensities and Zeeman 
frequency shifts are shown in figure I for "He (tt light 
polarisation). Shifts of transition frequencies are clearly 
visible on the projection onto the base plane, the usual 
Zeeman diagram (in which no line crossing occurs) . Large 
changes are also induced by the applied field on the transi- 
tion probabilities. In particular, the forbidden Ditransition 
(Y2 — > Z7) has a matrix element which linearly increases 

at low field (T2y''=3.3x |i?|), and approaches 1 in high field 
(curve labelled b in figure |^). The two other strong lines 
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Fig. 4. Level energies Es and Ep of the 2^S and 2^P states 
are computed as a function of the applied field B. They result 
from the fine-structure and Zeeman energy contributions for 
*He (dashed lines), with the additional hyperfine contributions 
for ^He (solid lines). The inserts are close-ups of the low-field 
regions. In the state (top), only the highest lying levels are 
shown for simplicity. 



at high field (2 curves labelled a superimposed in figure y) 
originate from the low-field Yi — > Z2 and Y3 — > Zg tran- 
sitions of D2 and Di lines. All other line intensities tend 
to at high field, in a way consistent with the sum rules 
of equation ^. 

The splittings of the Cg and Cg lines of '^He at mod- 
erate field are displayed in figure ^. Circularly polarised 
light, which is used in OP experiments to deposit angular 
momentum in the gas, is more efficiently absorbed for cr_ 
polarisation when a field is applied. This can tentatively 
be related to results of our OP experiments : at B=0.1 T, 
the most efficient pumping line is actually found to be 
Cg, a-. OP results at B—0.6 T reported in reference |^ 
also indicate that cr_ polarisation is more efficient, in spite 
of fluorescence measurements suggesting an imperfect po- 
larisation of the light. One must however note that OP 
efficiency does not only depend on light absorption, and 
that level structures and metastability exchange collisions 
(which will be discussed below) also play a key role. Let us 
finally mention that the forbidden Co transition becomes 




Fig. 5. 3D plot of the transition matrix elements for vr light 
polarisation in '*He, T^^ ^ (vr) , (vertical axis) and of transition 
energy difference e (with respect to Ci) as a function of B. To 
highlight the changes of the transition intensities, the symbol 
sizes (symbol areas) are proportional to T^j \ a and b are the 
strongest lines in high field (see text). The projection onto the 
base plane is the usual Zeeman diagram. 



allowed in an applied field, but that the transition prob- 
abilities remain very weak and vanish again at high field. 




e(T) 

Fig. 6. Plot of the Zeeman line shifts for the Cg and Cg lines 
of '^He. The symbol types refer to different light polarisation 
states. The symbol sizes (areas) are proportional to the corre- 
sponding transition matrix elements Tij. 
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2.1.3 Discussion 

To assess the consequences of the main approximation in 
our calculations, namely the restriction of the Hamilto- 
nian to the triplet configuration, we have made a sys- 
tematic comparison to the results of the full calculation 
including singlet-triplet mixing. 

In the case of ^He, exact energies are of course obtained 
for B=0, and level energy differences do not exceed 1 kHz 
for B=2 T. This very good agreement results from the 
fact that the lowest order contribution of the singlet ad- 
mixture in the state to the Zeeman effect is a third 
order perturbation ]30| , ^ . For similar reasons, computed 
transition probabilities are also almost identical, and our 
simple calculation is perfectly adequate for most purposes 
if the proper parameters of table 5| in the Appendix are 
used. 

For ^He we have compared the results of the full cal- 
culation involving all 24 sublevels in the 2P states as de- 
scribed by the Yale group and two forms of the sim- 
plified calculation restricted to the 18 sublevels of the 2^P 
state. The simpler form only contains the contact interac- 
tion term Apl ■ S in the hyperfine Hamiltonian of equa- 
tion ^. The more complete calculations makes use of the 
full hyperfine interaction, with the correction term Hj^Ys 
introducing two additional parameters, d and e. 

The simpler 1-parameter calculation is similar to the 
zero-field calculation of reference (with updated values 
of all energy parameters), and provides a limited accuracy 
as shown in figure 0. Errors on the computed energies of 
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Fig. 7. Errors in calculated energies for the 18 sublevels of the 
^He 2^P state resulting from the sole use of the contact term 
Apl ■ S in the hyperfine Hamiltonian are plotted for different 
magnetic field intensities. Large changes of the errors occur 
close to level anticrossings (e.g. mF = l/2 B15 and B17 levels 
close to 1.18 T, see figure §). 



the levels of several tens of MHz cannot be significantly re- 
duced by adjusting Ap to a more appropriate value. This 
is a direct evidence that the correction term H^Ys unam- 
biguously modifies the energy level diagram, with clear 



field-dependent signatures at moderate field and around 
1 T. Errors on the transition matrix elements up to a few 
lO^'^ also result from neclecting H^Ys- Using this simple 1- 
parameter form of the hyperfine Hamiltonian should thus 
be avoided when a better accuracy is desired. 

The 3-parameter calculation has also been compared 
to the full calculation, and a slight adjustment of param- 
eters Ap, d and e with respect to the corresponding pa- 
rameters C, D/2 and E/5 oi the full calculation (which 
also includes the off-diagonal parameter C |3^), now pro- 
vides a more satisfactory agreement for the level energies, 
shown in figure H. Transitions matrix elements Ty are also 
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Fig. 8. Errors in calculated energies for the 18 sublevels of the 
^He 2^P state resulting from the use of the full 3-parameter 
hyperfine Hamiltonian are plotted for different magnetic field 
intensities. As in figure ^, large changes of the errors occur 
close to the B15 - B17 level anticrossing around 1.18 T. 



obtained with good accuracy (LSxlO""* r.m.s. difference 
on average for all transitions and fields, with a maximum 
difference of 5x 10""*). The values of the 3 parameters Ap, 
d and e differ by -0.019%, 3.4% and 4.6% from the corre- 
sponding parameters in the full calculation. This adjust- 
ment may effectively represent part of the off-diagonal ef- 
fects in the singlet admixture, and allows to decrease the 
energy errors, especially below 1 T (by a factor 2-3). 

The other important approximation in our calculations 
is related to the simplifications in the Zeeman Hamiltonian 
written in equations ^ and || The first neglected term is 
the non-diagonal contribu tion proportional to the small 
parameter in equation The non-zero elements of 
this additional contribution are of order fisBgx/b, i.e. 1- 
2 MHz/T in frequency units. The other neglected term 
results from the quadratic Zeeman effect , which 

introduces matrix elements of order {^bB) /i?oo (-R00 — 
3.29 X 10^ GHz stands for the Rydberg constant in fre- 
quency units), i.e. 60 kHz/T^. For most applications, and 
in particular for the experimental situations considered in 
the following, corrections introduced by these additional 
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terms in the Zeeman Hamiltonian would indeed be nee 
gible, not larger for instance than the errors in figure p 



2.2 Metastability exchange 

So far we have only addressed the efi^ects of an applied 
magnetic field on the level energies and structures of the 
2^S and 2"^? states of helium, and the consequences for the 
1083 nm optical transition. In this section we now consider 
the effect of the applied field on the angular momentum 
transfer between atoms during the so-called metastability 
exchange (ME) collisions, which are spin-exchange binary 
collisions between an atom in the ground state and one 
in the metastable state. They play an important role 
in all experimental situations where the atoms in the 2'^S 
state are a minority species in a plasma (e.g. with a num- 
ber density 10^^-10^^ cm~'^ compared to 2.6x10"^^ for the 
total density of atoms at 1 mbar in typical OP conditions) . 
These collisions may be neglected only for experiments on 
atomic beams or trapped atoms for which radiative forces 
are used to separate the metastable atoms from the much 
denser gas of atoms in the ground state. 

To describe the statistical properties of a mixture of 
ground state and 2^S state atoms we shall use a standard 
density operator formalism and extend the treatment in- 
troduced by Partridge and Series for ME collisions, 
later improved and used for an OP model to 

the case of an arbitrary magnetic field. Our firstgoal is 
to show that the "spin temperature" concept |47|^ ,[7t re- 
mains valid, and that the population distribution in the 
2^S state is fully determined (in the absence of OP and of 
relaxation) by the nuclear polarisation M of the ground 
state. This is the res ult on which the optical detection 
method of section 3.3 relies. The second goal is to provide 



a formalism allowing to study the consequences of hyper- 
fine decoupling on spin transfer during ME collisions. 



2.2.1 Derivation of rate equations 

We shall consider situations in which no resonance is driven 
between the two magnetic sublevels of the ground state of 
^He (the ground state of ^He has no structure) nor be- 
tween sublevels of the 2'^S state. We can thus a priori 
assume that the density operators pg (in the ground state 
of ^He), p3 (in the 2^S state of ^He) and p4 (in the 2^S 
state of *He) are statistical mixtures of eigenstates of the 
Hamiltonian. The corresponding density matrices are thus 
diagonal and contain only populations (no coherences). 
With obvious notations : 



Pg 



l + M 
2 

^6 



1 - M 



where at and yi are relative populations (^ai—'^yi—l 



(11) 
(12) 
(13) 



An important feature of ME in helium is that in prac- 
tice no depolarisation occurs during the collisions, due to 
the fact that all involved angular momenta are spins [^6| . 
Three kinds of ME collisions occur in an isotopic mixture, 
depending on the nature of the colliding atoms (^He or 
''He) and on their state (ground state: '^'^He, 2'^S state: 
3.4He*). 

1. Following a ME collision between ^He atoms: 
3He + 3He* ^ ^He* + ^Re 



Pa 



P3 



P3 



the nuclear and electronic angular momenta are re- 
combined in such a way that the density operators just 
after collision p'^ and pg are given by: 



P'g ^^'^ P3 

P'z = Pg® Tr„ p3 



(14) 
(15) 



where Trg, Tr„ are trace operators over the electronic 
and nuclear variables respectively | |45[ | . 
2. Following the collision between a ground state '*He 
atom and a 2^S state ^He atom: 



^He 



3He* ^ ^ 3jje 

P3 P'L p'g 

the density operators just after collision are: 

P3=TreP3 (16) 
P4=Tr„p3. (17) 



3. Following the collision between a 2'^S state ''He atom 
and a ground state '^He atom: 



3He 

Pg 



"He* 

P4 



3Hc* 

n'" 
P3 



^He 



the density operator of the outgoing state ^He 
atom is: 



Pz =pg'i<> Pi- 



rn 



The partial trace operations in equations ^ to 
do not introduce any coherence term in density matri- 
ces. In contrast, the tensor products in equations ^ and 
introduce several off-diagonal terms. These are driving 
terms which may lead to the development of coherences 
in the density operators, but we shall show in the fol- 
lowing that they can usually be neglected. Since is 
a good quantum number, off-diagonal terms in the ten- 
sor products only appear between 2 states of equal mp. 
With the state names defined in the Appendix (table p^ , 
these off-diagonal terms are proportional to the operators 
and jA*^) (A'^l (for mF=l/2), \A''_) {AL\ and 
(for mi?=-l/2). These four operators contain 



A'^ 



AL) {A'l_ 

a similar sin6'±cos6'± factor, where 6+ and 6- are field- 



depe ndent level mixing parameters (see equations A14 to 
[A17 in the Appendix). They have a time evolution char- 
acterised by a fast precessing phase, exp{iAE±t/h) , de- 
pending on the frequency splittings AE±/h of the eigen- 
states of equal mp. 
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The time evolution of the density operators due to ME 
collisions is obtained from a detailed balance of the de- 
parture and arrival processes for the 3 kinds of collisions. 
This is similar to the method introduced by Dupont-Roc 
et al. |45| to derive a rate equation for pg and in pure 
^He at low magnetic field. It is extended here to isotopic 
mixtures, and the effect of an arbitrary magnetic field is 
discussed. 

The contribution of ME collisions to a rate equation for 
the density operator p^ of the 2^S state '^He is obtained 
computing an ensemble average in the gas of the terms 
arising from collisions with a ^He atom (collision type 1, 
first line of the rate equation) and from collisions involving 
different isotopes (types 2 and 3, second line): 



( 4-"-3P3 ) = {nzN^av^z {-ps + p^)) 



+ (aw34 {-n-iNips + n4iV3/53")> (19) 

in which the brackets ( ) correspond to the ensemble av- 
erage, iVa 4 and 713 4 are the ground state and 2'^S state 
atom number densities of ^He and "^He, a is the ME cross 
section, W33 and W34 are the relative velocities of collid- 
ing atoms. The ensemble averaging procedure must take 
into account the time evolution between collisions. For 
non-resonant phenomena, one is interested in steady-state 
situations, or in using rate equations to compute slow evo- 
lutions of the density operators (compared to the ME col- 
lision rate l/r). For the constant or slowly varying pop- 
ulations (y03, diagonal terms in pg and p'^'), the ensemble 
averages simply introduce the usual thermally averaged 
quantities 0^33^. For each time-dependent off-diagonal 
coherence, the average of the fast precessing phase involves 
a weighting factor exp(— t/r) to account for the distribu- 
tion of times i elapsed after a collision: 

1 f°° 

(exp {iAE±t/h)) = - / exp {iAE±i/h) 
''' Jo 

= {l-iAE±T/h)-^. (20) 

The ME collision rate l/r is proportional to the helium 
pressure (e.g. 3.75x10^ s~^/mbar in pure ^He at room 
temperature |4^). The splittings AE±/h are equal to the 
hyperfine splitting (6.74 GHz) in zero field, and t heir varr 



ation with B can be derived from equations AlO and All 



AE+/h increases with B while AE-/h initially decreases 
slightly (down to 6.35 GHz at 0.08 T), and both splittings 
increase linearly with B at high field (with a slope of or- 
der 27.7 GHz/T). Under such conditions, AE±T/h>1000 
in a gas at 1 mbar. All coherences can thus be safely ne- 
glected in equation [l^ which, when restricted to its diag- 
onal terms, can be written as: 



dt 



ME 



nsNsavss (-ps + ^ U^pg (g) Tr„ psU^ 



crw34 (^-n^Nips + UiNs ^ n,pg p4^i) (21) 



in which we note ilj =|Ai 
state Aj . 



(Ail the projector on the eigen- 



The main differences with the result derived in |45| are 
the additional term (second line) in equation |2j, which re- 
flects the effect of "^He— ^He collisions, and the replacement 
of the projectors on the F substates by the more general 
projectors Ui on the eigenstates in the applied field B. 
Considering the field dependence of the frequency split- 
tings AE±/h, the condition on the ME collision rate l/r 
is only slightly more stringent at 0.08 T that for B—Q, 
but can be considerably relaxed at high magnetic field. 
The linear increase of AE± and the 1 /B decrease of sin 9± 
(see figure ^l]) in the common factor of all coherences pro- 
vide a B^ increase of the acceptable collision rate, and 
hence of the operating pressure, for which coherences can 
be neglected in all density matrices, and equations |ll] to 

are valid. 

The other rate equations are directly derived from equa- 
tions |l|, [| and [l7 . Since there is no tensor product, hence 
no coherence source term in the relevant arrival terms, ex- 
plicit ensemble averages are not required and the contri- 
bution of ME is: 



dt 



ME 



n^Nsavss {~Pg + Trg p^) 
+ CTwil {-n^NsPg + n^N^ Tie pz) (22) 



— n4P4 ) = -UiNsavs^Pi + n^Niav^a Tr„ /93. (23) 
at J 



A trace operation performed on equations ^ and p3| shows 
that 



dt 



Hi 



ME 



dt 



"3 



(24) 



ME 



as expected since the three ME processes preserve the to- 
tal number of atoms in the 2'^S excited state. 



2.2.2 Steady-state solutions 

Total rate equations can be obtained by adding relaxation 
and OP terms to the ME terms of equations ^ and ^ 
They can be used to compute steady-state density ma- 
trices for the 2'^S state, while the nuclear polarisation M 
(and hence Pg) may have a very slow rate of change (com- 
pared to the ME collision rate l/r). In these steady-state 
situations, the contribution of ME collisions to the total 
rate equations for the 2^S state just compensates the OP 
and relaxation contributions. Since the latter are traceless 
terms (OP and relaxation only operate population trans- 
fers between sublevels of a given isotope), by taking the 
trace of equation |2^ or |2^ one finds that the 2^S state 
and ground state number densities have the same isotopic 
ratio R: 



R = Hi/ns = Ni/Na. 



(25) 



This results from having implicitly assumed in section 2.2.1 
that all kinds of collisions have the same ME cross sec- 
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tion|^ In fact destruction and excitation of of the 
state in the gas may differently affect '^He and ^He atoms, 
e.g. due to the different diffusion times to the cell walls. 
However in practice this has little effect compared to the 
very frequent ME collisions, which impose the condition 
of equation ^ for number densities. 

To obtain simpler expressions for the rate equations, 
we introduce the parameters: 



1/T, 



(26) 
(27) 
(28) 



The rates l/Tg and 1/re correspond to the total ME col- 
lision rate for a "^He atom in the ground state and in the 
2^S state respectively; in steady state, equation |2| can be 
used to show that Tg/Tg = n^/N^ = n^/Ni. Since the ther- 
mal velocity distributions simply scale with the reduced 
masses of colliding atoms, the value of the dimensionless 
parameter ^ can be evaluated from the energy dependence 
of tJ. From reference one can estimate /i ~ 1.07 at room 
temperature. Rewriting equation |2q as: 



1/re = [Nz + Ni)mj^ 



1 + R/n 
1 + R ' 



(29) 



the ME rate in an isotopic mixture thus differs from that 
in pure "^He by at most 7 % (at large R) for a given total 
density. 

Using these notations and the isotopic ratio R (equa- 
tion ^5|), the contributions of ME to rate equations be- 
come: 



dtP' 



ME 



1 1 



ME 



ME 



Tre ps) 

n , (-P4 + Tr„/93) 



(30) 
(31) 



X iiipg ® Tr„ p3 + Rpg ® p4) iTj . (32) 



In a way similar to that of reference we can trans- 
form equations ^ to ^ into an equivalent set of coupled 
rate equations for the nuclear polarisation M and for the 



^ The assumption that ME cross sections have the same value 
(7 for any isotopic combination of colliding atoms is valid at 
room temperature. This would be untrue at low temperature 
due to a 10 K isotopic energy difference in the 20 eV excitation 
energy of the 2^S state atoms (the energy is lower for '^He). 



relative populations ai and yi of the 2 S states and Y,;: 

6 \ 



dt 



di^' 



dt 



ME 



ME 



ME 




1 1 



TeR 
1 



-a., - 



k.ak 



Eg; 



(33) 
(34) 



1 



R + p 

3 

MFl)ak + RY.{E:, 



k=l 



fc=l 



(35) 



The values of the M-independent but B-dependent ma- 
trices L, G^, E^, and F^ are provided in the 
Appendix. Equation p3 ^ d irectly results from computing 
Tr„ Pglz using equation |30|. The rate equations on the rela- 
tive populations are obtained by computing (Y^ | dp4^/dt JY^) 
using equation ^ and (Ai| dp^/dt \Ai) using equation |2[ 
The linear A/-dependence in equation |3^ directly results 
from that of pg (equation |ll| ) . 



2.2.3 Spin temperature distributions 

Anderson et al. have proposed that, under some condi- 
tions including fast spin exchange, the relative populations 
of sublevels should follow a Boltzmann-like distribution in 
angular momentum. This was explicitly verified for pure 
■^He in low field 0| , and will now be shown for isotopic mix- 
tures and arbitrary magnetic field. In situations where OP 
and relaxation processes have negligible effect on popula- 
tions, the steady-state density operators are easily derived 
from the rate equations M to ^: 



Pg = Tre P3 
P4 = Tr„ p3 



(36) 
(37) 

(38) 



Simply assuming that the populations ai in the 2^S state 
of ^He only depend on the mp values in th e state s A^, 
one can directly check using equations A14 to A17 that 
both Tre Ps (and hence equations ^ andpsj) and Tr„ p^ 
do not depend on B. Noting e'' the ratio (l+A//)/(l-M) 
of the populations in the ground state (1//3 plays the role 
of a spin temperature), one derives from equation ^ that 
the ratios of populations in the 2^S state of '^He are field- 
independent, and given by e^'^™" . These populations thus 
have the same distribution at all magnetic fields: 



/(e3/5/2 + 2e^/2 + 2e-^/^ + e'^f^l^). (39) 



Using equations M and A14 to A17, the populations in 
the state of ^e are in turn found to obey a similar 
distribution: 



1 



(40) 
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2.2.4 OP effects on populations 

The general problem of computing all the atomic popu- 
lations in arbitrary conditions is not considered in this 
work. It has been addressed in the low field limit using 
specific models in pure ■^Hc |^] and in isotopic mixtures 
|p(]| . Here we extend it to arbitrary magnetic fields, but 
only consider the particular situation where there is no 
ground state nuclear polarisation (M=0). It can be met 
in absorption spectroscopy experiments, in which a weak 
probing light beam may induce a deviation from the uni- 
form population distribution imposed by M=0 (infinite 
spin temperature). We shall assume in the following that 
the OP process results from a depopulation mechanism in 
which population changes are created only by excitation 
of atoms from selected sublevels and not by sponta- 
neous emission from the state (where populations are 
randomised by fast relaxation processes). In ^He, this de- 
population mechanism has been checked to dominate for 
pressures above '--^l mbar |^j7|. 

Under such conditions, and further assuming that Zee- 
man splittings are large enough for a single sublevel in 
the 2'^S state to be pumped (Ap or Yp, depending on the 
probed isotope), the OP contribution to the time evolu- 
tion of populations has a very simple form. Given the ME 
contribution of equation ^ or ^5|, the total rate equation 
for the pumped isotope becomes: 



d 








d 




dt"' = 





ME 



ME 



<5y(«,p) 



5A{i-p) 



(41) 
(42) 



where I/t^ is the pumping rate (equation ^ or ||) and 
^a(*,p) = 1 (resp. 5Y{i,p) — '^) if the level (resp. Y^) is 
the pumped level, otherwise. In steady state, the relative 
populations can be obtained from the kernel of the 9x9 
matrix representation of the set of rate equations ^ and 
H, or m and |l|. 

When a ■'He line is probed, the ''He relative popula- 
tions are fully imposed by ME processes only (equa- 
tion Q), and can be subsituted in equation |4^. The rela- 
tive populations can then be obtained as the kernel of a 
6x6 matrix which depends on the pumped level index p, 
on the reduced pumping rate Te/Tp and on the magnetic 
field S, but not on the isotopic ratio R (as discussed at 
the end of the Appendix). Results showing the decrease 
of the pumped population a^, and hence of the absorp- 
tion of the incident light, are displayed in figure ^ as a 
function of i?. A strong absorption decrease is found at 
high B, due to larger population changes induces by OP 
when ME less efficiently transfers the angular momentum 
to the ground state (let us recall that we assume M—Q, 
so that the ground state is a reservoir of infinite spin tem- 
perature). Population changes are found to be weaker for 
the two states A*^ and A'i, which both have |0 : ±) as the 
main component in high field (curves b and b' in figure 

I). 

When a He line is probed, the substitution and re- 
duction performed above cannot be made. The full set of 
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Fig. 9. OP on a ''He level : the computed decrease of the 
relative population a-p of the pumped level of '^He below its 
equilibrium value (1/6) is plotted as a function of B for two 
values of the reduced pumping rate Te/rp : 10"'^ (solid lines) 
and 10"'^ (dashed lines). For the 2 populations ap=a\ or a'l, 
the results fall on the same curve (labelled b or b' depending 
on T^/Tp). For the 4 other populations, a significantly larger 
decrease is found (groups of curves a and a'). The insert is an 
expanded view of the low- field region (with a linear B-scale). 



populations turns out to vary with the isotopic ratio R. 
However, the OP effects are found to depend mostly on 



the product RTf,/Tp^ as shown in figure |10| in the case of 
the state Y3 (7715=!). This influence of the isotopic com- 
position extends over a wide field range, which justifies the 
use of isotopic mixtures to reduce the bias resulting from 
OP effects in our systematic quantitative measurements 
of "^He line intensities. Figure |ll| displays a comparison 
of computed OP effects for the 3 sublevels of the ''He 2^S 
states. Results for the populations yi and are quite sim- 
ilar, except at large R. Weaker OP effects are found for 
the population 2/2, for which 7715=0. This feature is similar 
to the weaker OP effects found in '^He for the two states 
which have m5=0, m/=±l/2 as main angular momenta 
components in high field (see figure ||). 



3 Experimental 

In this part we present experimental results of laser ab- 
sorption experiments probing the fine, hyperfine and Zee- 
man splittings of the 1083 nm transition in ^He and ''He. 
Accurate microwave measurements [p7|j33j34| a nd high pre- 
cision laser spectroscopy experiments |28H29| , ^1| , ^ are usu- 
ally performed on helium atomic beams. In contrast, we 
have made unsophisticated absorption measurements us- 
ing a single free-running laser diode on helium gas in cells 
at room temperature, under the usual operating condi- 
tions for nuclear polarisation of '^He by OP. An important 
consequence of these conditions is the Doppler broadening 
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6(1) 

Fig. 10. OP on a *He level : the decrease of the *He population 
y3 below its equilibrium value (1/3) is plotted as a function 
of B. The population has been computed for two values of 
the reduced pumping rate Te/rp : (solid lines) and 10~^ 

(dashed lines), and different values of R (ranging from 100 to 
10""^). The population decrease mostly depends on the product 
Rre/rp (10"\ 10"^, ...10"'^ for curves a, b, ... e). 




e(T) 

Fig. 11. OP on a ''He level : the computed decrease of ''He 
populations yi, j/2 and j/3 due to OP effects is plotted as a 
function of B. The reduced pumping rate Te/rp ^lO"'^ is that 
of the solid lines in figure lid. Four values of the isotopic ratio 
R are used (Rre/rp =10" to 10"** for curves a to d). 



of all absorption lines. Each transition frequency LUij /2t: 
(in the rest frames of the atoms) is replaced by a Gaussian 
distribution of width A : 



(43) 



depending on the atomic mass Mat- At room tempera- 
ture (T=300 K), the Doppler widths for ^He and -^He 
are Zig^Lig GHz and zi4=1.03 GHz (the fuU widths at 
half maximum, FWHM, given by 2Z\\/Tn2, are respec- 



tively 1.98 and 1.72 GHz). The narrow Lorentzian res- 
onant absorption profiles (with widths r' 12) in equations 
H and [s] are thus strongly modified, and broader Voigt 
profiles are obtained for isolated lines in standard absorp- 
tion experiments. For F' ^ Z\, these are almost Gaus- 
sian profiles of width A. In the following, we first report 
on Doppler-free absorption measurements from which a 
good determination of frequency splittings is obtained, 
then on integrated absorption intensities in different mag- 
netic fields. We finally demonstrate that absorption mea- 
surements may provide a convenient determination of the 
nuclear polarisation M . 



3.1 Line positions: saturated absorption measurements 

3.1.1 Experimental setup 



The experiment arrangement is sketched in figure 12. The 
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^abry Perot interferometer 
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Fig. 12. Main elements of a saturated absorption experiment. 
A magnetic field B is applied either along the pump and probe 
beams, or perpendicular to the beams. Polarising elements 
(Pol.) consist of linear polarisers and/or quarter-wave plates 
depending on the desired polarisations of the pump and probe 
beams. 



helium gas samples are enclosed in sealed cylindrical Pyrex 
glass cells, 5 cm in diameter and 5 cm in length. Cells are 
filled after a careful cleaning procedure: bakeout at 700 K 
under high vacuum for several days, followed by 100 W RF 
(27 MHz) or microwave (2.4 GHz) discharges in helium 
with frequent gas changes until only helium lines are de- 
tected in the plasma fluorescence. Most results presented 
here have been recorded in cells filled with 0.53 mbar of 
■^He or helium mixture (25% ■^He, 75% ''He). Different gas 
samples (from 0.2 to 30 mbar, from 10% to 100% ^He) 
have also been used to study the effect of pressure and 
isotopic composition on the observed signals. A weak RF 
discharge (<1 W at 3 MHz) is used to populate the 2^S 
state in the cell during absorption measurements. The RF 
excitation is obtained using a pair of external electrodes. 
Aligning the RF electric field with B provides a higher 
density of 2'^S states and better OP results in high mag- 
netic field. 

Most data have been acquired using a specially de- 
signed air-core resistive magnet (100 mm bore diameter. 
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with transverse optical access for 20 mm beams). In spite 
of a reduced footprint (30x20 cm^, axis height 15 cm) 
which conveniently permits instaUation on an optical ta- 
ble, this magnet provides a fair homogeneity. The com- 
puted relative inhomogeneity is <10~^ over a 5 cm long 
cylindrical volume 1 cm in diameter (typical volume probed 
by the light beams in spectroscopy measurements), and 
<3xl0^^ over the total cell volume. This usually induces 
low enough magnetic relaxation of nuclear polarisation in 
OP experiments. The field-to-current ratio was calibrated 
using optically detected NMR resonance of '^He. Experi- 
ments for B up to 0.12 T were performed with this coil. 
Several saturated absorption experiments have been re- 
peated and extended up to 0.22 T using a standard elec- 
tromagnet in the Institute of Physics in Krakow. 

The laser source is a 50 mW laser diode (model 6702- 
Hl formerly manufactured by Spectra Diode Laborato- 
ries). Its output is coUimated into a quasi-parallel beam 
(typically 2x6 mm in size) using an anti- reflection coated 
lens (/==8 mm). A wedge-shaped plate is used to spht 
the beam into a main pump beam (92% of the total in- 
tensity), a reference beam sent through a confocal Fabry 
Perot interferometer, and a probe beam (figure 12). A 
small aperture limiting the probe beam diameter is used to 
only probe atoms lying in the central region of the pump 
beam. A small angle is set between the counterpropagating 
beams. This angle and sufficient optical isolation from the 
interferometer are required to avoid any feedback of light 
onto the laser diode, so that the emitted laser frequency 
only depends on internal parameters of the laser. The po- 
larisations of pump and probe beams are adjusted using 
combinations of polarising cubes, 1/2- wave, 1/4- wave re- 
tarding plates. For all the measurements of the present 
work, the same polarisation was used for pump and probe 
beams. However valuable information on coUisional pro- 
cesses can be obtained using different polarisation and/or 
different frequencies for the two beams. 

The probe beam absorption is measured using a modu- 
lation technique. The RF discharge intensity is modulated 
at a frequency /rf low enough for the density of the ab- 
sorbing atoms 2^S to synchronously vary (/rf ^^lOOHz). 
The signal from the photodiode monitoring the transmit- 
ted probe beam is analysed using a lock-in amplifier. The 
amplitude of the probe modulation thus measured, and 
the average value of the transmitted probe intensity are 
both sampled (at 20 Hz), digitised, and stored. Absorp- 
tion spectra are obtained from the ratios of the probe 
modulation by the probe average intensity. This proce- 
dure strongly reduces effects of laser intensity changes and 
of optical thickness of the gas on the measured absorp- 
tions n. 

With its DBR (Distributed Bragg Reflector) technol- 
ogy, this laser diode combines a good monochromaticity 
and easy frequency tuning. It could be further narrowed 
and stabilised using external selective elements, and thus 
used for high-resolution spectroscopy of the 1083 nm line 
of helium pl| . Here we have simply used a free running 
diode, with a linewidth dominated by the Schalow-Townes 
broadening factor of order 2-3 MHz FWHM H. The laser 



frequency depends on temperature and current with typi- 
cal sensitivities 20 GHz/K and 0.5 GHz/mA respectively. 
To reduce temperature drifts, the diode is enclosed in a 
temperature regulated 1-dm'^ shield. Using a custom made 
diode temperature and current controller, we obtain an 
overall temperature stability better than a mK, and a 
current stability of order a fiA. Frequency drifts over sev- 
eral minutes and low frequency jitter are thus limited to 
10 MHz at most. 

The laser frequency is adjusted, locked or swept by 
temperature control using the on-chip Peltier cell and tem- 
perature sensor. The small deviation from linearity of the 
frequency response was systematically measured during 
the frequency sweeps by recording the peaks transmitted 
by the 150 MHz FSR interferometer. Figure [l^ displays an 
example of such a measurement: for a control voltage V (in 
Volts) the frequency offset (in GHz) is 52.5 x (1+0.048^)1/. 
The linear term in this conversion factor is usually taken 
from more accurate zero-field line splitting measurements 
(see below), but this small non-linear correction is used 
to improve the determination of frequency scales for ex- 
tended sweeps. 
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Fig. 13. Insert ; fraction of recording of the light intensity 
transmitted by the Fabry Perot interferometer during a tem- 
perature controlled laser frequency scan (1 V of control voltage 
approximately induces changes of 2.5 K in laser temperature 
and 50 GHz in emitted frequency). Main plot : residues of a 
linear fit (open triangles) or parabolic fit (solid circles) of all 
transmitted peak positions. 



3.1.2 Zero-field measurements 

Results presented in this section have in fact been ob- 
tained in the earth field in which Zeeman energy changes 
are of order 1 MHz or less and thus too small to affect the 
results given the accuracy of our measurements. A small 
permanent magnet is placed near the cell to strongly re- 
duce the field homogeneity and thus prevent any nuclear 
polarisation to build up when the pump beam is absorbed. 
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An example of frequency sweep over the whole '^He ab- 
sorption spectrum is given in figure |l^ ('^He pressure: 0.53 
mbar). The frequency scale is determined using the ac- 
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Fig. 14. Zero-field saturated absorption spectrum of '^He. The 
upper traces (a) represent the probe absorption signals when 
the pump beam is applied (solid line) or blocked (dotted line) . 
Traces b and c (an expanded section of trace b) are the differ- 
ence of probe absorptions with and without the pump beam, 
with a x4 vertical scale factor. The sections of the spectrum 
comprising the lines C2 to C5 provide a much larger signal and 
are displayed with a reduced amplitude. Frequency sweep rate 
is 160 MHz/s, a value ensuring sufficient signal quality and 
manageable laser frequency drifts. Transition lines are labelled 
Ci to C9, crossover resonances between Ci and Cj are labelled 
Xi-j- The vertical bars on the axes represent the computed 
line positions and intensities. 



curately known value 3As/2=6.7397 GHz of the Cs-Cg 
spUtting. The saturated absorption spectrum (solid line in 
figure [l4|a) combines the usual broad absorption lines and 
several narrow dips. The dips reveal a reduced absorption 
of the probe by atoms interacting with the pump beam. 
These are atoms with negligible velocity along the beam 
direction (Doppler-free resonances) or atoms with a veloc- 
ity such that the Doppler shift is just half of the splitting of 
two transitions from or to a common level (crossover reso- 
nances). An absorption signal recorded without the pump 
beam (dotted line in figure |lja) only shows the broad fea- 
tures, which can be fit by a Doppler Gaussian profile for 
isolated lines (e.g. Cs or Cg). ft is used to obtain the nar- 
row resonances in figure Hjb by signal substraction. The 
width S of the narrow resonances is 50 MHz FWHM in 
the conditions of figure nj. ft can be attributed to several 



combined broadening processes, including pump satura- 
tion effects and signal acquisition filtering. 

Broader line features, the "pedestals" on which nar- 
row resonances stand, result from the combined effects of 
collisions in the gas. Fully elastic velocity changing col- 
lisions tend to impose the same population distribution 
in all velocity classes, and thus to identically decrease 
the probe absorption over the whole velocity profile. In 
contrast, ME collisions involving a '^He atom in a gas 
where the nuclear polarisation is M~0 induce a net loss 
of angular momentum^ (a uniform population distribu- 
tion tends to be imposed by the infinite spin tempera- 



ture, see section 2.2.3). In steady-state, the un-pumped 
velocity classes will thus acquire only part of the pop- 
ulation changes imposed in the pumped velocity class. 
The relative amplitude of the resulting broad pedestal 
results from the competition of coUisional transfer and 
collisional loss of orientation. One consequence is that 
it is almost pressure- independent, as was experimentally 
checked. Another consequence is that it strongly depends 
on ^He concentration in isotopic mixtures, as illustrated 
in figure 15l The much larger fractional amplitude of the 
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Fig. 15. Zero-field saturated absorption spectrum of the Do 
line in a 25% '^He - 75% *He mixture (same total pressure, 
0.53 mbar, as for figure The lower trace also displays the 
squared absorption amplitude (open circles), with an arbitrary 
weight adjusted to fit the saturated absorption pedestal (see 
text). 



pedestal results from the weaker depolarising effect of the 
ME collisions: only 25% of the exchange collisions occur 
in this case with a depolarised '^He atom and thus con- 



The orientation is not totally lost at each ME collision since 
the electronic part of the angular momentum is conserved, and 
recoupled after collision. ME collisions thus contribute both to 
orientation transfer between velocity classes and to orientation 
loss. The relative importance of these contributions depends 
on the angular momentum loss, which is reduced in high field 
by hyperfine decoupling. A specific feature of ME collisions is 
that they involve small impact parameters and large collision 
energies due to a centrifugal energy barrier p"8|j2^ , and thus 
usually produce large velocity changes. 
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tribute to bleach the pump-induced population changes, 
while all collisions induce velocity changes. In figure ^ 
is also displayed the squared absorption profile, which 
closely matches the shape of the pedestal. This is a general 
feature also observed in ^He spectra, which results from 
the linear response to the rather low pump power. Assum- 
ing that the relative population changes are proportional 
to the absorbed pump power, frequency detuning reduces 
by the same Doppler profile factor both these population 
changes and the probed total population. The resulting 
squared Doppler profile is still Gaussian, but with reduced 
width Z\/^/2. 

From this detailed analysis of the signal shapes we in- 
fer two main results. First, very little systematic frequency 
shift is expected to result from the distortion of the nar- 
row Doppler-free resonances by the broad line pedestals of 
neighbouring transitions (no shift at all for isolated lines). 
It can be estimated to be of order (5^//\, i.e. of order 1 MHz 
and thus negligible in these experiments. Second, popula- 
tion changes in the whole velocity profile remain limited 
in pure '^He or in '^He-rich isotopic mixtures. They amount 
for instance to a 7 % change for the Cg transition when the 
pump beam is applied, as illustrated in figure H (differ- 
ence between dotted and solid line) . They are significantly 
larger in a helium mixture, e.g 30 % for the conditions of 
figure |l^. These OP effects on the populations would be 
reduced for a reduced pump power, and the direct effect 
of the attenuated (-^25) probe itself can be expected to be 
correspondingly reduced, and hence negligible under simi- 
lar conditions. However, this would be untrue in pure *He 
in which strong population relaxation is difficult to im- 
pose, and precise absorption or saturated absorption ex- 
periments are less conveniently performed. A strong mag- 
netic field would also enhance the pertur bing e ffects of the 
probe beam, as was disc ussed in section 2.2.4 and will be 
demonstrated in section 3.2.4. 



A series of identical sweeps including the C1-C7 lines, 
similar to that in figure ^4|c, was performed to check for 
reproducibility. The frequency scales were determined us- 
ing the C2-C6 splittings (these transitions connect the two 
hyperfine levels of the 2^S state to the same excited level, 
see figure ^, and are thus separated by 6.7397 GHz just as 
Cs-Cg). To directly probe hyperfine level splittings in the 
state we consider the line separations C2-C4 and C3- 
C5, which directly measure the level energy differences be- 
tween the lower pairs of 2^P levels (see the zero field level 
diagram in figure^). We also consider the separation of the 
crossover resonances X2-4-X3-5, actually computed from 
differences of the four line positions, which measures an 
independent energy difference between these levels. Com- 
parison to the values computed using the 3-pa rameter hy- 
perfine term in the Hamiltonian (section ^.1.3| ) is shown in 
figure |l6|. For the two simple splittings, data scatter is ±15 
MHz about the average. It is smaller for the third com- 
bined splitting, which involves averaging pairs of line posi- 
tion measurements and thus a statistically reduced uncer- 
tainty. The three averages have a 6 MHz r.m.s. difference 
with the computed values, which is consistent with statis- 
tical uncertainty. This analysis shows that our unsophisti- 
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Fig. 16. Difference between measured and computed hyperfine 
splittings in the 2''P state of '^He in zero field for selected pairs 
of individual or combined line splittings (see text). Each closed 
symbol corresponds to a different experiment. The large open 
symbols are the statistical averages. 



cated spectroscopic measurements are quite accurate and 
have very little systematic error in the frequency deter- 
minations (at most a few MHz over several GHz). It also 
justifies the use of the correction term beyond the contact 
interaction term in the hyperfine contribution (equation ||) 
since our experimental errors are lower than the effect of 
this correction term (figure |^). 



3.1.3 Field effect on line positions 

The most obvious effect of an applied magnetic field is 
to considerably increase the number of observed narrow 
lines (Doppler-free lines and intercombination lines) in the 
saturated absor ptio n recordings. For simplicity, we only 
present in figure|l7| selected results for the set of '^He lines 
originating at low field from Cg and Cg. These are the 
most efficient optical transitions for OP of '^He, and thus 
are of special interest. Furthermore, the analyses of tran- 
sitio n in tensity measurements described in sections 3.2.2 
and 3.3 request an accurate knowledge of all line positions 
for these transitions. To reduce the number of observed 
lines and thus facilitate data processing, the same circu- 
lar polarisation (cr-|_ or (t_) was used for the pump and 
the probe beam, and two measurements were thus per- 
formed for each field value. For the higher field data (iron 
yoke magnet, open symbols in figure ^^a) the use of an 
imperfect circular polarisation made it possible to usually 
observe all intense lines in each measurement. Zero field 
scans were performed between measurements to evaluate 
the slow drift of the laser frequency and provide reference 
measurements for the Zeeman shifts. 

For the measurements in the air-core magnet (up to 
0.125 T), the agreement between computed and measured 
line positions e/h is fair, but affected by fluctuating fre- 
quency offsets. The r.m.s. differences are quite large (50 
MHz and 110 MHz for a+ and ct_ probes respectively), 
but the largest contribution to these differences can be 
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Fig. 17. a : Plot of frequency shifts on the He Cg and Cg 
lines in an applied magnetic field measured by a saturated ab- 
sorption technique. Solid lines : computed Zeeman shifts (same 
as in figure ^ . Up and down triangles : data for cr+ and cr_ po- 
larisations of light beams. Solid and open symbols : measure- 
ment using the air-core and iron yoke magnets respectively. 
Crosses : result of corrections applied the open symbol data 
set (see text), b : Differences between experimental and com- 
puted values for the measurements in the iron yoke magnet 
are plotted as a function of B. Open symbols : raw data (no 
correction) ; filled symbols : data with field corrections only ; 
crosses : data with field and frequency corrections (see text). 



attributed to a global frequency offset from scan to scan. 
Estimating the offset of each scan from the average of the 3 
recorded Une positions, the remaining differences are sig- 
nificantly reduced (9 MHz r.m.s. for ct+ and ct_). The 
accuracy of these line position measurements is thus com- 
parable to that of the zero-field measurements (figure p^ . 

For the measurements in the iron yoke electromagnet 
(up to 0.25 T), a similar frequency offset adjustment is 
not sufficient to significantly reduce the differences be- 
tween measured and computed line positions (figure [T^b) . 
For instance, they amount to 130 MHz r.m.s. for the a+ 
measurements (open up triangles). Moreover, splittings 
within a given scan consistently differ from the computed 
ones. Assuming that the magnetic field is not exactly pro- 
portional to the applied current due to the presence of a 
soft iron yoke in the magnet, we use the largest measured 
splitting on the Cg line (the transitions originating from 
the sublevels mi?=±3/2) to determine the actual value of 
the magnetic field during each scan. Solid symbols in fig- 
ure p7[ b are computed differences for these actual values 
of the magnetic field. Crosses (figures |l7| a and [l7|b) corre- 
spond to fully corrected data, allowing also for frequency 



offset adjusments. The final remaining differences (22 MHz 
r.m.s.) are still larger than in the air-core magnet, yet this 
agreement is sufficient to support our field correction and 
offset adjustment procedures in view of the line intensity 
measurements described in the following section. 



3.2 Line intensities 

3.2.1 Experimental setup 

The experimental arrangement is either the one sketched 
in figure 12 , in which the pump beam is blocked and only a 
weak probe beam (usually 0.1-1 mW/cm^) is transmitted 



through the cell, or a simplified version of the setup when 
no saturated absorption measurement is performed. 

The magnetic field is obtained by different means de- 
pending on its intensity. Data at moderate field (up to 
0.22 T) data have been acquir ed using the resistive mag- 
nets described in section 3.1.1. Higher field measurements 
have been performed in the bore and fringe field of a 1.5 T 
MRI superconducting magnet |5^. With the laser source 
and all the electronics remaining in a low-field region sev- 
eral meters away from the magnet bore, cells were succes- 
sively placed at five locations on the magnet axis with field 
values ranging from B=6 mT to 1.5 T. The 1.5 T value 
(in the bore) is accurately known from the routinely mea- 
sured ^H NMR frequency (63.830 MHz, hence 1.4992 T). 
The lower field intensities, 6 mT and 0.4 T, have been mea- 
sured using a Hall probe with a nominal accuracy of 1%. 
The intermediate values, 0.95 T and 1.33 T, are deduced 
from the measured Zeeman splittings with a similar accu- 
racy of 1%. Due to the steep magnetic field decrease near 
the edge of the magnet bore, large gradients cause a sig- 
nificant inhomogeneous broadening of the abso rption lines 
in some situations discussed in section 3.2.3. For these 



high field measurements, the Fabry Perot interferometer 
was not implemented and no accurate on-site check of the 
laser frequency scale was performed. Instead, the usual 
corrections determined du ring t he saturated absorption 
measurements (see section 3.1.1 and figure |l^) are used 
for the analysis of the recorded absorption spectra. 

A systematic study of the effect of experimental con- 
ditions such as gas pressure, RF discharge intensity, and 
probe beam power and intensity, has been performed both 
in low field (the earth field) and at 0.08 T. The main char- 
acteristic features of the recorded spectra have been found 
to be quite insensitive to these parameters. Absorption 
spectra have been analysed in detail to accurately deduce 
the relative populations of different sublevels (from the ra- 
tios of line intensities) and the average atomic density of 
atoms in the 2^S state in the probe beam (from its absorp- 
tion). A careful analysis of lineshapes and linewidths re- 
veals three different systematic effects which are discussed 
in the next three sections. 



3.2.2 Effect of imperfect ^He isotopic purity 

Evidence of a systematic effect on absorption profiles was 
observed in experiments with nominally pure '^He gas. A 
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typical recording of the Cs and Cg absorption lines for 
B=0 and residuals to different Gaussian fits are shown 
in figure A straightforward fit by two independent 
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Fig. 18. Zero-field absorption measurement on the Cs and 
Cg lines in nominally pure ^He at 0.53 mbar. a : absorption 
signal (the vertical scale is 20 times less sensitive than for the 
other traces), b : residue plot of independent Gaussian fits to 
the two lines, c : residue plot of the same fit as in plot b for 
Cg and the corresponding computed Cg profile, d : computed 
D2-D1 profile with a -^200 amplitude scaling corresponding 
to the measured traces of ^He (see text), e : residue plot of 
independent Gaussian fits to the two lines after correction for 
the *He contribution. 



Gaussian profiles to the recorded lines (trace a) system- 
atically provides a 1.5-2 % larger width and an unsatis- 
factory fit for the Cg transition (trace b, residue plot). 
Moreover, the ratio of fitted line strengths (areas) is 7.8% 
larger than the value computed using table ^ in the Ap- 
pendix (1.374 instead of 1.274). Conversely, fitting on the 
Cg component and assigning the expected lincwidth, am- 
plitude and frequency shift to the Cg component provides 
the residue plot of trace c, which suggests the existence 
of an additional contribution to absorption. All this is ex- 
plained by the presence of traces of ^He, which are in- 
dependently observed by absorption measurements on the 
isolated Dq line. Indeed this small amount of ^He (isotopic 
ratio R=0.2 to 0.5 %) affects the absorption measurement 
performed on '^He due to the intense lines which lie within 
one Doppler width of the Cg line (trace d: computed ab- 
sorption lines of ^He for R=0.2 %). The reduction of the 
^He gas purity is attributed to the cell cleaning process, 
which involves several strong discharges performed in '*He 
prior to final filling with pure '^He. Reversible penetration 
of helium into the glass walls is believed to be responsi- 
ble for the subsequent presence of a small (and surpris- 
ingly unsteady) proportion of "^He in the gas. When the 
appropriate correction is introduced to substract the con- 
tribution of "'He lines from the absorption spectra, the fit 
by two independent Gaussian profiles is significantly im- 
proved (trace e) and the linewidths are found to be equal 



within ±0.2 %. For the recording of figure [l^a one obtains 
a Doppler width Z\=1.198 GHz, corresponding to a tem- 
perature of 304 K in the gas (consistent with the ambient 
temperature and the power deposition of the RF discharge 
in the plasma). 

The ratio of the line amplitudes after this correction 
is closer to the ratio of computed transition intensities, 
still 3.5±0.2% larger in this example. A series of measure- 
ments performed in a row provides very close results, with 
a scatter consistent with the quoted statistical uncertainty 
of each analysis. In contrast, a measurement performed in 
the same cell after several days may lead to a different 
proportion of ^He and a different ratio of transitions in- 
tensities (for instance i?=0.35% and a ratio of intensities 
2.9±0.1% smaller than the computed one). 

Consequences of the presence of a small proportion of 
^He on the lineshapes and intensities of recorded spec- 
tra are indeed also observed in an applied magnetic field, 
as shown in figure |l^ for i3=0.111 T. The recorded ab- 
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Fig. 19. Absorption measurements on the Cg and Cg lines 
performed in a field _B=0.111 T in nominally pure ^He at 0.53 
mbar for circular polarisations (upper graph : cr+ , lower graph 
: (T_). Traces a correspond to the raw absorption signals (solid 
lines) and computed profiles (dotted lines, see text). Traces 
b are differences (residue plots) between absorption signals, 
corrected for light polarisation defects, and computed profiles. 
Traces c are the computed profiles for pure *He, scaled down 
in order to correspond to an isotopic ratio _R=0.4%. Expanded 
vertical scales (x5, right hand side axes) are used for traces b 
and c. 



sorption signals are similar to the computed ones (dot- 



18 



E. Courtade et al.: The 1083 nm Helium Transition in High Magnetic Field 



ted lines), but two significant differences appear. First, 
a small line is visible for instance in the upper graph 
around e/h—30 GHz, where no (t+ transition should be 
observed. This is due to an imperfect circular polarisation 
of the probe beam in the cell, resulting in part from stress- 
induced circular dichroism in the cell windows (measured 
to be of order 0.2%). The amount of light with the wrong 
polarisation is determined for each recording, and each sig- 
nal is corrected by substracting the appropriate fraction 
of the signal with the other circular polarisation (0.79% 
and 0.18% corrections for a+ and a- recordings in fig- 
ure |l^)0. Second, the corrected signals may then tenta- 
tively be fit by the sum of three Gaussian profiles centred 
on the frequencies independently determined in a satu- 
rated absorption experiment (see figure But, as is the 
case in zero field, this procedure does not provide satisfac- 
tory results for the amplitudes and linewidths of the split 
Cg lines (corresponding to the transitions originating from 
any of the lowest levels of the 2''S state, Ai to A4, to the 
highest levels of the 2^P state, B17 and Big). In figure 
traces b display the differences between the corrected sig- 
nals and the computed sums of three Gaussian profiles 
with a common linewidth and with amplitudes propor- 
tional to the nominal transition probabilities in this field. 
Two global amplitudes and the linewidth (Z\=1.221 GHz) 
are adjusted to fit the two Cg transitions, and the resulting 
differences are best accounted for assuming a proportion 
i?=0.4% of '*He in the gas (traces c in figure ^9^. In order 
to better reproduce the experimental recordings, the am- 
plitudes of the split Cg lines are finally allowed to vary. 
Fits with excellent reduced E^^id no visible feature in 
the residues are obtained and the obtained ratios of in- 
tensities are given in table ^ for this typical recording . 
The same kind of discrepancy as in zero field is obtained 
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Table 1. Computed line intensities Tij for _B=0.1111 T and 
their ratios are compared to the experimental line amplitude 
ratios of figure piJ. 



for the strong components of the transitions; it happens 
to be much worse for the weakest line in some recordings, 
including that of figure ^ for reasons which are not un- 
derstood. 

In an attempt to test whether the remaining line in- 
tensity disagreements with the computed values may re- 



* Imperfect light polarisation has no effect in very low field 
where Zeeman energy splittings are negligible and transition 
intensities of each component (Ci to Cg and Dq to D2) do not 
depend on light polarisation. 



suit from imperfect corrections for the presence of "^He, or 
from difficulties due to fitting several overlapping lines, we 
have also probed the well resolved Dq line of ''He in a gas 
mixture in the 0-0.12 T field range. Using a circularly po- 
larised probe beam, the two transitions connecting levels 
Yi and Y3 to Zg (cr^ and cr_ polarisations re spect ively) 
are split in frequency by 56 GHz/T (see section ^.1.2| ). The 
ratio of the corresponding line intensities is computed to 
be approximately given in this field range by: 



^it) w 1 + 1.88B + 2.02^2, 



-(4) 
.9 



(44) 

where B is the field intensity in Tesla. The results of a 
series of absorption measurements for a transverse probe 
beam (with linear perpendicular polarisation) are displayed 
in figure EOl The absorption signals are fit by the sum of 
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Fig. 20. Absorption measurements on the Do line of ''He per- 
formed in fields up to 5=0.12 T with a lin+ polarised probe 
beam in a ^He-^He mixture at 0.53 mbar (isotopic ratio _R=3). 
a : Recordings of absorption signals are given for two values of 
the field. The lower plots display the results of Gaussian fits 
to the signals, b : Doppler width A4, c : ratios of measured to 



computed relative line amplitudes (A^^/A^-) / (73*gV7i\^ 



,9 



three Gaussian profiles centred on the transition frequen- 
cies (measured in a saturated absorption experiment). This 
allows for a small amount of tt polarisation in the probe 
beam, which experimentally arises from imperfect align- 
ment of the direction of the polarisation with respect to 
B. In spite of this drawback, the transverse probe beam 
configuration is preferred to a longitudinal one, for which 
polarisation defects may directly affect the ratio of line 
amplitudes. 

The common Hnewidth Zi4=1.067 ±0.006 GHz (fig- 
ure corresponding to a temperature of 321 K) is found 
to be independent of the field. The ratio Acr+/Ao-^ of the 
measured line amplitudes for the cr+ and a- polarisation 
components is found to increase with B as expected from 
equation p4|, with however a slightly lower value (^--^3%, 



figure 20 



To summarise the results of this section, the differ- 
ence of measured line intensities in our experiments with 
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respect to the computed ones is usually small (of order 
±3%) after correcting for the presence of some ^He in our 
^He cells and for the observed light polarisation imperfec- 
tions. This difference does not arise from signal-to-noise 
limitations and its origin is not known. This may set a 
practical limit to the accuracy with which an absolute 
measurement of population ratios can be made using this 
experimental technique. Similar discrepancies are also ob- 
served at higher fields, situations for which the additional 
systematic effects described in the following sections must 
first be discussed. 



3.2.3 Effects of magnetic field gradients 



As previously mentioned (section |3.2.1 ), absorption ex- 
periments above 0.25 T have been performed in the fringe 
field of a 1.5 T magnet. The exact field map of this magnet 
was not measured, but the field decreases from 1.33 T to 
0.95 T in only 25 cm, so that field gradients along the axis 
may exceed 15 mT/cm at some locations. In contrast, field 
variations in transverse planes are much smaller in the 
vicinity of the field axis. As a result, a negligible broad- 
ening of absorption lines is induced by the field gradient 
for perpendicular probe beams (cr or tt polarisation for 
a linear polarisation perpendicular or parallel to B). For 
probe beams propagating along the field axis {a polarisa- 
tion for a linear polarisation, or tT_ polarisation for a 
circular polarisation), the range of fields SB applied to the 
probed atoms is the product of the field gradient by the 
cell length. This field spread SB induces spreads Seij/h 
of optical transition frequencies which are proportional to 
SB and to the derivatives dtij/dB. Each spread Scij/h 
thus depends on the probed transition as illustrated in 
figure P for 5=0.95 T. 



quency spread functions, computed from the density dis- 
tribution of metastable atoms, the field gradient at the cell 
location and the derivatives dtij/dB. Assuming for sim- 
plicity a uniform density of atoms in the 2^S state, com- 
puted broadened absorption profiles only depend on the 
dimensionless ratio Scij/hA of the transition frequency 
spread to the Doppler width. Figures |2^ a and |2^ display 
examples of computed profiles (absorption as a function of 
this reduced detuning). For increasing frequency spreads. 
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Fig. 22. a : Computed absorption (half-)profiles in a mag- 
netic field gradient resulting from the convolution of a Gaussian 
Doppler profile (width A) and of a uniform spread of the res- 
onance frequency (over a width Seij/h). b : The same profiles 
are plotted in semi-logarithmic scale as a function of the square 
of the detuning. A Gaussian profile appears as a straight line, 
c : Apparent linewidths Aapp and amplitudes Aapp of these 
computed profiles are plotted as a function of the square of 
the reduced resonance frequency spreads. The solid symbols 
correspond to the profiles displayed in a and b. 
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Fig. 21. Computed field-derivatives of the transitions energies 
eij/h are plotted for all ^He and ''He circular polarisation tran- 
sitions as a function of eij/h for _B=0.95 T. Inhomogeneous line 
broadenings proportional to these values are induced by a field 
gradient. 

The resulting absorption profiles are given by the con- 
volution of the Doppler velocity distribution and of the fre- 



profiles become broader and their maximum amplitude de- 
creases (the area is strictly independent of Seij/h A). De- 
viation from a true Gaussian profile is more clearly seen in 
figure p2| b as a deviation from a straight line. To test the 
procedure used in the analysis of experimental data, one 
may fit a Gaussian function to these computed profiles. 
The apparent linewidths and amplitudes thus obtained, 
Aapp and Aapp^ are plotted in figure [2^c. For small fre- 
quency spreads, Aapp and 1/Aapp are almost proportional 
to Sefj. However, significant deviations appear at large 
Scij, and the product AappAapp is not strictly constant 
(it increases by 2 % for Seij/hA ~3). This means that 
whenever frequency spreads cannot be avoided (e.g. using 
a transverse probe beam), fitting procedures to Gaussian 
profiles introduce systematic errors. Computation of line 
areas would provide the true line intensities, but this pro- 
cedure cannot be used to separately obtain the intensities 
of ill-resolved transitions. 

Figure ^ displays part of an absorption measurement 
performed at 5=0.95 T with a longitudinal cr-polarised 
probe beam. It includes the two split Dq transitions con- 
necting levels Yi and Y3 to Zg : the a+ ^He line with 
highest e/h (101 GHz), and (t_ ^He line with highest 
e/h (47 GHz), respectively. The recorded absorption sig- 
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Fig. 23. Part of an absorption spectrum recorded at 5=0.95 T 
in a '^He-*He mixture at 2.13 mbar (isotopic ratio R=3) using 
a longitudinal probe beam (open symbols). The solid curve is 
the spectrum computed for this R, assuming thermal Doppler 
widths for the different lines. Its amplitude is scaled to repro- 
duce that of the 47 GHz line. Experimental broadening of the 
absorption lines is directly found on the 101 GHz line, and 
suggested by the residual absorption between the main lines 
(non-zero signal at 50 and 55 GHz), which is 2-3 times larger 
than computed. The insert is an expanded portion of the spec- 
trum, in which several '^He and ''He lines overlap (the computed 
absorption of pure ''He is displayed as a dotted line). Experi- 
mental broadening of the absorption lines also clearly appears 
for this part of the spectrum. 




((e-e(4)..)//7)2 (GHz2) 

Fig. 24. Semi-logarithmic plot of absorption data recorded at 
5=0.95 T in a '^He-^He mixture at 2.13 mbar (isotopic ratio 
R—3) . They are presented for 3 different *He lines as a function 

of the square of the frequency mistuning — e'^'^ /h. Open 
symbols : data of figure ^ longitudinal probe beam. Solid 
symbols : data for a transverse probe beam. The straight lines 
represent Gaussian absorption profiles, and the values of their 
slopes are used to extract the apparent linewidths Aapp- The 
smallest linewidths, consistent with the 'He Doppler width, are 
obtained for a transverse probe (Z\4=1.13 and 1.105 GHz for 
the 47 and 72 GHz lines). Larger apparent linewidths are ob- 
tained for a longitudinal probe {Aapp=l-27, 1.31 and 1.42 GHz 
for the 47, 72 and 101 GHz lines). 



nal (symbols) roughly correspond to the profile computed 
assuming the nominal Doppler widths A and neglecting 
frequency spreads (solid line). Clear amplitude differences 
are however observed, in particular for the most intense 
line. They can be attributed to OP e ffects, and will be 
discussed in detail in the next section (3.2.4). Moreover, 
significant linewidth differences appear, for instance on the 
isolated most shifted ^He line (fitted apparent linewidth 
Aapp=lA2 GHz, 30% larger than A4). Most of the other 
lines actually result from the overlap of several atomic 
transitions, and a straightforward Gaussian fit is not rele- 
vant. For example, the computed total absorption around 
e/h=72 GHz (insert in figure |2^, solid line) results from 
the two *He transitions at 66.7 and 72.1 GHz (dotted line) 
and the two ^He transitions at 65.2 and 71.1 GHz. 

To better analyse the effect of the probe beam direc- 
tion on the shape and widths of different lines, we plot 
in figure ^ the absorption data using the same method 
as in figure 223. The open symbols correspond to three 
of the lines in figure |2^, with a probe beam along the di- 
rection field, and the solid symbols to lines from a similar 
recording with a perpendicular probe beam. Apart from 
the asymmetries on the 71 GHz lines which result from 
the additional '^He line on the low-frequency side of the 
absorption profile (which is thus not retained in the anal- 
ysis), satisfactory linear fits reveal Gaussian-like profiles 
and provide the values of the corresponding linewidths. 
The experiments performed with a transverse probe pro- 
vide data with the steepest slopes and hence the smallest 



linewidths, while those using a longitudinal probe lead to 
larger apparent linewidths as expected . Using the convo- 
lution results of figure p2| c, the frequency spreads Seij/h 
for these lines are computed to be 1.4, 1.65 and 1.85 GHz 
respectively. From the derivatives dcij/dB for these transi- 
tions (figure |l]) , one computes that a field spread SB=52 mT 
would produce frequency spreads of 1.08, 1.65 and 1.85 
GHz, corresponding to the observed broadenings for the 
72 and 101 GHz lines. The additional apparent broaden- 
ing for the 47 GHz experimental line may be attributed 
to O P effe cts for this rather intense line (see the next sec- 
tion 3.2.4| ). For the estimated value of the field gradient 
(15 mT/cm), this 52 mT field spread would be obtained 
for an effective cell length of 3.5 cm. This is quite satis- 
factory since the density of atoms in the metastable 2^S 
state vanishes at the cell walls, so that most of the probed 
atoms lie within a distance shorter than the actual gap 
between the end windows of our 5 cm-long cells (external 
length) . 

To summarise the results obtained in this section, the 
expected effects of strong field gradients on absorption 
lines have indeed been observed. They introduce a sys- 
tematic broadening and amplitude decrease of absorption 
signals which can be significant (up to 30 % changes in 
our experiments), and must be considered in quantitative 
analyses. For isolated lines, area measurements are not af- 
fected by the transition frequency spreads. However, most 
lines usually overlap in a recorded absorption spectrum 
(especially in an isotopic mixture), and the quantitative 
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analysis of such spectra is in practice too inaccurate if the 
different hnewidths are not known, ft is thus important 
to only use a transverse probe beam configuration when- 
ever accurate absorption measurements are needed and a 
field gradient is present. The same argument would also 
apply to saturated absorption measurements, for which 
even limited frequency spreads (e.g. of order 100 MHz) 
would strongly decrease the signal amplitudes. 



3.2.4 Effects of optical pumping 



line 


d (Do, a-) 


Do (B=0) 


a (Do, (7+) 


Tig 


0.7722 


0.3333 


0.06459 


weight (calc.) 


19.31 % 


8.333 % 


1.615 % 


weight (exp.) 


16.6 % 


7.9 % 


1.8 % 



Table 2. Computed line intensities Tig of the Do line for B=0 
and 1.5 T. The weights in the spectra (fraction of the line 
area over the total area) are computed for an isotopic ratio 
R=3 corresponding to the experimental conditions. The last 
line displays the experimental weights for the measurements 
shown in figure E3. 



A comparison of absorption recordings at B—Q and 1.5 T 
is presented in figure UM. As expected in such a high mag- 
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Fig. 25. Recordings of complete absorption spectra at 5=0 
and 1.5 T for circular polarisations in a ^He-'*He mixture at 
2.13 mbar (isotopic ratio _R=3). The lines labelled a to e are 
especially considered in the analysis of OP effects. The *He 
transitions with the highest energy (involving the Zg level) 
considered in figure ^ for moderate fields are labelled a and 
d at 1.5 T, and correspond to the Do line at zero field. The 
assignment of the lines at zero field can be made using figure |^. 



netic field, the absorption spectra are deeply modified and 
strongly depend on the probe light polarisation. Line shifts 
of order of the fine structure splittings, and large changes 
in the intensities of the absorption lines are observed, for 
instance on the well-resolvcd transitions to the excited Zg 
level (labelled a, Dq and d in figure which we shall 
generically call Dq transitions). For these lines, the com- 
puted transition intensities and the corresponding weights 
in the total spectra are given in table |[ It is not possible to 



directly compare absorption signals at B=0 and 1.5 T be- 
cause the plasma intensity and distribution in the cell are 
strongly modified in high field. However, according to the 
sum rules of equation S, the average density of atoms in 
the 2^S state along the probe beam is deduced by integrat- 
ing the absorption signal over the whole spectrum. This is 
used to scale the absorption amplitudes in figure and 
have a meaningful comparison of line weights in table |^. 
While there is a satisfactory agreement with expectations 
for all line weights at B=0, the ratio of the weight of the 
strong (Dp, cr_) absorption hue d to that of the weak (Dp, 
cr-(-) line a is substantially smaller than computed. 

The measurements displayed in figure ^ have been 
performed with a rather intense probe beam, so that OP 
effects may play a significant role. As for all experiments 
performed in low field, we can assume that OP creates no 
nuclear polarisation M for the ground state '^He atoms. 
It is not possible to enfor ce rap id nuclear relaxation as is 
done in low field (section 3.1.2) by imposing a strong rel- 
ative field inhomogeneity at 1.5 T. However the angular 
momentum deposited in the cell by the absorption of a 
fraction of the probe beam could only lead to a very slow 
growth rate of M for our experimental condition^ As a 
result, only a very small nuclear polarisation M can be 
reached in a frequency sweep during the time spent on an 
absorption line (~5 s). Negligible changes of the Zeeman 
populations and hence of absorption signal intensities are 
induced by this process. In contrast, the OP e ffects in the 
2^S metastable state described in section 2.2.4 may have a 



strong influence on the distribution of the Zeeman popu- 
lations, and thus significantly decrease the absorption (see 
figure 10). 

To discuss and analyse in detail these OP effects at 
1.5 T, we show in figure 26 (symbols) parts of the exper- 
imental spectra extracted from figure together with 
different computed profiles which are discussed below. It 
is important to note that in spite of the double-pass config- 
uration in the cell used to enhance absorption signals, no 
narrow saturated absorption dips appear in the data. This 
is different from what is observed at low field when a coun- 



^ It can be estimated from the number of '^He nuclei in the 
cell (10^^ for a partial pressure of 0.53 mbar) and the absorbed 
photon flux from the light beam (e.g. 4x10^'* photons/s for 
0.1 mW, corresponding to a 10 % absorption of a 1 mW probe 
beam): the growth rate of Al is of order 0.05 %/s. This result is 
consistent with purnping rates measured for optimal OP con- 
ditions at low field [p5[). 
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Fig. 26. Expanded plot of selected absorption profiles from ex- 
perimental recordings at 1.5 T (open symbols). The labelling 
of the lines is that of figure |2^. The vertical scale has been 
reduced (xO.3) for the strongest line d (the actual ratio of the 
maximum absorptions of lines d and a is 7.0). Dotted lines : 
computed Gaussian profiles for these transitions. Solid lines : 
computed profiles taking into account OP effects with the pa- 
rameters given in the text (hardly distinguishable from the 
dotted Gaussian profiles for lines a and b). 



terpropagating pump beam induces population changes 
in a given velocity class, and when velocity changing colli- 
sions only partly redistribute these changes over the whole 
Doppler profile (see section 3.1.2 and figures |lj - |l5|). 
We shall thus assume that in the conditions of figure pS] 
no correlation is introduced by the probe beam between 
atomic velocity and Zeeman population changes, i.e. that 
any out-of-equilibrium orientation produced by OP due to 
the probe is uniformly redistributed among all atoms in 
the Doppler profile. The pumping time Tp introduced in 
equations ^ and ^ can thus be computed from the pho- 
ton absorption rates defined in equations ^ and ^ using an 
averaging procedure similar to that of 0| . The averaging is 
performed over the atomic velocity projection v along the 
probe beam, which introduces the Doppler shift luv/c for 



the atomic transition angular frequencies uj=uj[^^ and 
in equations H and ^. For "^He transitions, one obtains^: 



1 



(r72)2e-(''/«)'dv 



(45) 



where the mean velocity v = 2'ircAy,lujij is related to 
the Doppler width defined in equation |4^. For a single- 
frequency laser and a small enough damping rate of the 
optical transition coherence [F' /2 ^ 2ttA^), the integral 



^ This differs by a factor l/\/2 from the (incorrect) corre- 
sponding result in reference |7|. 



is easily computed: 

1 0ra/ 



Fij lias 



(46) 



A similar result is obtained for He transitions by substi- 
tution of Z\4, T^'p and w^j^^ in equation The numerical 
value of the first factor appearing in the pumping rate 
expression is for ^He: 



nipUjAA 



= 4.269 X 10^ s"V(WAn ), 



(47) 



and it is 15 % smaller for '^He due to the larger value of 
the Doppler width. 

For a given laser intensity, the pumping rate is thus 
expected to be proportional to the transition intensity Tij 
and, for each absorption line, to decrease with increas- 
ing detuning from transition frequen cies o f atoms at rest. 
Relying on the OP model of section 2.2.4, Zeeman popu- 
lations can be numerically computed as a function of the 
pumping rate. The absorption signals, which are propor- 
tional to the probed populations, to and to the Gaus- 
sian factor in equation can then in principle be com- 
puted using a single global scaling parameter (the number 
density of atoms in the 2^S state). 

In our experiments, the laser intensity lias of the probe 
beam at the cell location could not be precisely measured, 
and instead was deduced from the comparison of the mea- 
sured ratio of the maximum absorptions of the cr_ and cr^ 
Do lines, labelled d and a in figures Eq and 
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to the com- 
puted ratio which would be obtained in the absence of 
OP effects, i.e. for yi=?/3=l/3 (y3T3^g/?/iTi^9=7.0 instead 
of T3,9/ri,9=11.96). The value (re/r^) /T^ =3.28x IQ-^ is 
numerically obtained from this comparison of OP effects 
on these lines at resonance (the two values for Tij are given 
in table 0). The exchange rate l/rg is computed using 
equation and the value of the exchange rate l/r^ in a 
pure ■^He gas at the operating pressure(3.75x 10^ s~^/mbar 



i9 



l/re = l/r°(l + i?(/i-l)). 



(48) 



For the conditions of the experiment (isotopic ratio i?=3, 
total pressure 2.13 mbar) re = 1.04xl0~^ s, hence l/{TpTij) 
=3.15x10"* s and Iias=7A W/m^ or 0.74 mW/cm^, which 
is quite consistent with its independently estimated value. 
For this value of lias , OP effects are computed for all lines 
as a function of the detuning from resonance, and the 
resulting absorption profiles precisely agree with the ex- 
perimental data (figure |2^, solid lines, compared to the 
Doppler profiles, dotted lines). This agreement over the 
whole profiles is better seen in figure E^. There is a slight 
systematic deviation for all lines at large detuning, where 
absorption is larger than computed. This may be due to a 
temperature slightly above room temperature (as is sys- 
tematically observed in other experiments), and to the 
difference between the true atomic (Voigt) profile and 
a Gaussian one. The agreement is otherwise very good 
and fully supports the coUisional redistribution assump- 
tion used to derive equation hn. These measurements can 
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Furthermore, the simple OP model of section 2.2.4 can- 
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Fig. 27. Semi-logarithmic plot of absorption data (open sym- 
bols) recorded at 5=1.5 T in a ''He-''He mixture at 2.13 mbar 
(isotopic ratio 7i=3). They are presented for 3 different ^He 
lines (Do, 3 light polarisations) as a function of the square of 
the frequency detuning e/h— e^^^ /h. The dotted straight lines 
represent gaussian absorption profiles, the solid lines are the 
computed profiles with OP effects (see text). 



also be viewed as indirect checks of the variation of the 
population of the pumped sublevel with the pumping rate, 
such as plotted in figures ^ and ^ up to a 45 % popula- 
tion reduction for 7/3. The influence of lias was indeed also 
directly measured by using a linear a polarisation instead 
of a circular one, thus reducing the effective pump power 
by a factor 2 on isolated lines. All measured OP effects 
are reduced according to calculations, with for instance a 
29 % maximum population reduction for on the (T_ Do 
line. 

OP effects for other transitions also agree with expec- 
tations. For instance, probing the population 1/2 with a 
TT polarisation (figure induces relatively small OP ef- 
fects in spite of a significant transition intensity. This is 
consistent with the reduced sensit ivity of the Y2 state to 
OP which was noticed in section 2.2.4| (figure |ri|). Most 
^He transitions are difficult to analyse due to the proxim- 
ity of several overlapping lines, but the lines c and e in 
figure ^ can be considered to evaluate OP effects. Line 
c is an isolated line probing the A5 state, so weak that 
no significant OP effect is expected (a 0.2 % maximum 
population reduction for 05), and observed. In contrast, 
line e consists of 4 lines (two with e/h ~-16.1 GHz, two 
with e/h ^--^-14. 9 GHz), each with a transition intensity 
~1. Experimentally, the ratio of the maximum absorp- 
tions of lines e and c is 40.2, 4 % lower than expected 
from transition intensities, thus suggesting a 4 % relative 
decrease of the probed populations. However these transi- 
tions concern 4 different sublevels in the 2^S state, this is 
an average decrease for the populations in these sublevels. 



not be directly applied since it assumes that a single level 
is pumped. Still, if the A5 state alone was pumped with 
the probe intensity used in the experiment, a computed 
3.2 % maximum population reduction would be obtained 
for 05, in reasonable agreement with the small observed 
OP effects. 

Let us finally mention that OP effects have been ob- 
served to have all the expected features at lower magnetic 
field. For instance, for the same laser intensity lias as at 
1.5 T, the maximum population reduction for y^, on the 
(Do, cr-) line is 38 % at 1.3 T, 26 % at 0.95 T and 8 % at 
0.4 T. The much weaker effects expected for other transi- 
tions could not clearly be observed, especially in the pres- 
ence of strong field gradients which also affect the appar- 
ent line intensities and shapes as discussed in the previous 
section. 

To summarise the results of this section, the computed 
effects of optical pumping in the 2^S state have been ob- 
served and quantitatively verified in some detail. This sup- 
ports the validity of the assumptions made on the effect 
of metastability exchange collisions on the Zeeman pop- 
ulations and of the simple OP model considered in sec- 
tion 2.2. It also confirms that orientation is fully redis- 



tributed by collisions over the whole velocity profile in 
high field. This can be understood if one considers that 
ME collisions usually strongly modify the velocities of the 
colliding atoms, but that due to hyperfine decoupling lit- 
tle angular momentum is transferred to the nuclear spin. 
An intuitive measure of this decoupling can be deduced 
from the high efficiency of OP processes at small pump- 
ing rate: a 45 % population decrease can be obtained at 
1.5 T for Te/Tp=l/400, which means that relaxation of 
electronic orientation typically requires 400 exchange colli- 
sions (or ~100 collisions with an unpolarised ground state 
•^He atom) . One consequence is that performing saturated 
absorption measurements similar to those of section 5.1 
to probe ^He transitions becomes increasingly difficult in 
fields above 0.5-1 T, and certainly requires the presence of 
a larger amount of '^He than in our experiments (an iso- 
topic ratio i? ^1). Finally, a practical result of this study 
of OP effects is that the attenuation level of a probe beam 
required to perform an unperturbed measurement of rel- 
ative populations is well determined. It strongly depends 
on the probed isotope and on the chosen transition, and 
on the gas composition and pressure, but a laser intensity 
lias typically 100 times lower (i.e. 10 /iW/cm^) than the 
highest one we have used must be chosen to have a 1 % 
accuracy on all population measurements up to 1.5 T. 



3.3 Optical measurement of nuclear polarisation 

In this section we address the measurement of the nuclear 
polarisation M of the ground state of '^He, such as can 
be efficiently prepared by performing OP of the 2'^S state, 
both in pure '^He and in isotopic mixtures. 

NMR is a valuable technique which can be used to 
directly measure the nuclear magnetisation resulting from 
the polarisation. It is however usually not well suited to 
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systematic studies of OP dynamics and efficiency: it has 
a reduced sensitivity for a low pressure gas, it induces 
polarisation losses at each measurement and thus cannot 
continuously measure M, and it is affected by the noise- 
generating plasma discharge used for OP. In addition, it 
requires a high absolute homogeneity of the magnetic field 
B, and a retuning or redesign of the NMR spectrometer 
at each change of the field value. 

In contrast, optical measurement techniques have a 
very high sensitivity to changes in the electronic states 
of the excited states of helium in the plasma, even in a 
low pressure gas. The standard optical detection technique 
|p|,p^ mentioned in the introduction, in which the circu- 
lar polarisation of a chosen helium spectral line emitted 
by the plasma is measured and the nuclear polarisation is 
inferred is simple to implement and provides an accurate 
non-destructive measurement of M. It relies on hyperfine 
coupling to transfer angular momentum from nuclear to 
electronic spins in the excited state which emits the mon- 
itored spectral line. The decoupling effect of an applied 
magnetic field unfortunately reduces the efficiency of this 
angular momentum transfer, which is also sensitive to de- 
polarising collisions. The usefulness of this technique is 
thus restricted to low fields (<10 mT), low gas pressures 
(<5 mbar) and limited '^He isotopic ratios {R <1) to avoid 
a significant sensitivity loss (-^2 for each of the quoted 
limits). Different optical methods, which rely on absorp- 
tion measurements on the 2^S - 2'^P transition, have been 
successfully used to quantitatively determine the nuclear 
polarisation of ^He P,|8|, |56| , ^ . It is a method of this kind, 
designed to operate for arbitrary magnetic fields and for 
all isotopic mixtures, which we describe and discuss in the 
rest of this article. 



3.3.1 Principle 

Optical absorption measurement methods provide infor- 
mation both on the total number density of atoms in the 
2'^S state and on the relative populations of the probed 
sublevels. In usual situations, the population distribution 
in the 2^S state is strongly coupled by ME collisions to 
that in the ground state, characterised by M. These pop- 
ulations would exactly be rul ed by a spin temperature 
distribution (see section 2.2.3 ) in the absence of OP or 
relaxation processes. For simplicity, we shall assume that 
these additional processes have a negligible effect com- 
pared to that of ME. The validi ty of t his assumption will 



be discussed in the final section 3.3.3 



When two absorption measurements directly probe two 
populations of atoms in the 2'^S state (for instance at low 
field when the line Cg or Dq is probed with cr+ and cr_ cir- 
cular polarisations), the derivation of M is a straightfor- 
ward procedure. When transitions simultaneously probe 
several sublevels (e.g. in low field with a polarisation on 
any line, with any polarisation on line Cg, etc.), the 
measurements of two independent combinations of pop- 
ulations can still be used to infer the nuclear polarisation 
M, but specific calculations are then required Here, 
we are mostly interested in measurements at high enough 



magnetic fields for the Zeeman shifts to remove all level 
degeneracies {B >50 mT). It is thus very easy to selec- 
tively probe a single level by an adequate choice of the 
transition frequency and of the polarisation of the light. 



With the spin temperature 1//3 defined in section 2.2.3 
the nuclear polarisation is given by: 



M - (e^ - 1) / (e^ + 1) 



(49) 



Assuming that the two measured absorption signals, S 
and 5', are simply proportional to the number densities of 
atoms in the probed sublevels (averaged in the volume of 
the probe beam), we can relate them to the total number 
density and to the relative populations. When ^He atoms 
are probed, one obtains: 



S — Kn^a and S' ~ K'n^a , 



(50) 



where only the total number density 713 and the relative 
populations a and a' of the two probed sublevels vary with 
M . The Af-independent coefficients K and K' depend on 
the probed transitions; they are proportional to the transi- 
tion intensities and thus depend on the chosen transitions 
and they vary with the magnetic field. The ratio S' / S 
only depends on the ratio of transition intensities and of 
relative populations (see equation B3) : 



S'/S ^ K'/Ke 



l3(m'p—mF) 



(51) 



where mp and m'p are the total angular momentum pro- 
jections in the two probed states. To eliminate systematic 
errors resulting from an inaccurate determination of the 
ratio K' / K, it is often convenient to also measure absorp- 
tion signals, 5*0 and S'q, in the unpolarised system {M—0, 
hence S'q/Sq — K'/K). The spin temperature can thus 
be directly deduced from ratios of measured absorption 
signals: 



5^ = {S'So/SS'o)'^^""'^-"'^^ 



(52) 



provided that m'p 7^ mp- The nuclear polarisation can 
then be obtained using equation The total number 
density 713 in the 2'^S state can in turn be computed from 
equations ^ ^ and In the following we only discuss 
its relative variation with M: 



n3{M) _ S e^^ + 2e^ + 2e-^ + e'^^ 
n3(0) ^ 6e^ ' 



(53) 



Similar results are indeed obtained when "^He atoms 
are probed to measure two relative populations, y and y', 
among the three sublevels in the 2''S state. In this case. 



the angular momentum projections in equations Bl 



^2| must be replaced by ms and m'g, and equation 53| is 
replaced by: 



and 



1 



ni{M) _ 

^4(0) ~ 3e'3™s 



(54) 
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3.3.2 Examples of operation 

As a demonstration of the performance of an absorption- 
based optical measurement of M in non-standard condi- 
tions, we report here on experimental measurements per- 
formed at 0.0925 T. This field value is high enough for 
the Zeeman splittings to be at least of the order of the 
Doppler line widths, so that generic problems arising from 
high- field operation are present. The main difference with 
low-field techniques is that the two absorption measure- 
ments required to record the signals S and S' must be 
performed using different probe laser frequencies (a few 
line crossings at some field values, e.g. around 0.16 T - see 
figure |l^ - may be used to avoid this constraint, but this 
will not be considered in the following). Before describing 
experimental schemes to obtain such dual-frequency mea- 
surements, we present examples of continuous frequency 
scans of the probe laser in steady-state OP situations. 

Figure |28| a displays recordings of absorption signals on 
the Cs and Cg '^He lines for a tt polarisation of the trans- 
verse probe. They are obtained with the sa me te chnique 
as for line intensity measurements (section |3.2.1 ), using 
a modulation (at a frequency /rf ^100 Hz) of the RF 
discharge excitation of the plasma and a phase-sensitive 
analysis at /rf of the transmitted probe intensity. The 
recorded spectra are quite similar to those of figure 
but four lines are obtained for this tt polarisation. A fit 
of the spectra by a sum of four Gaussian lines is used to 
obtain the relative line intensities given in table ^. As for 



e/h (GHz) 


30.61 


34.24 


38.60 


40.60 


level, mp 


Ag, + 


As, - 


A3, -f 


A2, - 


So/Tij 


1.003 


0.991 


0.989 


1.018 




0.890 


1.098 


0.884 


1.127 


S/So 


0.888 


1.109 


0.894 


1.108 


combination 




ae/a2 


as /as 


a-i/a'z 


M for So 


0.006 


-0.001 


-0.007 


-0.014 


M for S 


-0.105 


-0.108 


-0.117 


-0.121 


M for S/So 


-0.111 


-0.107 


-0.110 


-0.107 



Table 3. Table of normalised relative line intensities So/Tij 
and S/Tij obtained from fits of the two recordings of figure bqa. 
The nuclear polarisation M is computed from ratios of inten- 
sities for the four relevant combinations of lines {m'p 7^ mp). 
From the ratios S/ So (last line) , the nuclear polarisation is de- 
duced to be M=— 0.109±0.002. This negative polarisations is 
obtained due to the use of cr_ OP light. 



the absorption measurements discussed in section 3.2.2, 
several systematic effects actually limit the accuracy of 
the measurements of relative line intensities (e.g. imper- 
fect light polarisation or ■^He purity). A ±1.5 % devia- 
tion from the computed intensity ratios is observed for 
the M=0 recording in figure |l^ (first line in table ||). This 
accounts for the observed scatter in the determinations of 
M deduced from the ratios of signal amplitudes ( "A/ for 
S" values in table ^), and justifies the use of the absorp- 
tion measurements at null polarisation and of equation M 
to actually compute M (last line in table ||) . 





35 40 

e/h (GHz) 

Fig. 28. a : Experimental recordings of absorption signals ob- 
tained at 0.0925 T in pure '^He (pressure 8.5 mbar) for M=0 
(solid line) or in a steady-state OP situation (open symbols). 
The OP is performed using a laser diode (difi'erent from the one 
used to probe absorption) tuned on the maximum absorption 
of the ((J_, C9) transition. The incident pump power is 30 mW, 
and the transmitted beam is back reflected once through the 
cell. Lock-in detection is performed at the modulation fre- 
quency /rf of the RF discharge, b : Experimental recordings 
of signals obtained with a frequency-modulated probe beam 
and second-harmonic phase-sensitive detection (see text) in the 
same conditions. 



Replacing a continuous frequency scan by operation at 
two fixed frequencies requires accurate independent set- 
tings of the laser frequencies to probe the centres of the 
absorption lines. With a different modulation scheme, in 
which the laser frequency is modulated over a small frac- 
tion of the Doppler width at a frequency /probe ^1 kHz, a 
phase-sensitive analysis at /probe of the transmitted probe 
intensity provides an error signal proportional to the deriva- 
tive of the absorption signal. This error signal allows to 
lock the laser frequency on any absorption line. A second 
harmonic phase-sensitive analysis (at 2/probc) provides a 
signal proportional to the second derivative of the absorp- 
tion signal. This is displayed in figure p8| b for the same 
conditions as those of the absorption measurement in fig- 
ure 28 1. The four transitions appear to be better resolved, 
and the amplitude variations can similarly be used to de- 
duce the polarisation M. The respective advantages of the 
two modulation schemes will be discussed in the next sec- 
tion. 
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In isotopic mixtures the same transitions could still 
be probed, but the overlap of the intense "^He lines pre- 
cludes any reliable analysis of the Cg transition. Instead 
we present recordings of absorption spectra on the Dp lines 
in figure 29 for a transverse probe beam with equal weights 
on a and tt polarisations. While the same laser power is 




zlh (GHz) 

Fig. 29. Experimental recordings of absorption signals 
0.0925 T in an isotopic mixture (i?=3, total pressure 8.5 mbar) 
for M=0 (solid line) or in a steady-state OP situation (open 
symbols), a and tt polarisations are equally present in the 
transverse probe beam. The OP is performed in the same way 
as for the recordings in figure but on the (a_ , Do) transition. 
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Fig. 30. Nuclear polarisation M and relative changes in total 
metastable density na deduced from successive measurements 
of absorption signals by states A5 and A3 (_B=0.0925 T, '^He 
pressure 1.07 mbar). The OP is performed on the (cr+, Cg) 
transition with 150 mW of incident pumping power starting 
from an unpolarised gas (this sets the initial time t=Q). The 
pumping beam is blocked after 180 s and M exponentially 
decays back to (Ti=156 s). Here polarisation is positive due 
to the choice of a+ OP light. 



used for OP, a much larger nuclear polarisation (M=0.27) 
is deduced from this measurement in an isotopic mixture. 
This higher efficiency is similar to that systematically ob- 
served in low field OP experiments and is related to 
the more efficient absorption of the pumping light by the 
((T_,Do) fine. 

Frequency scans over 20 GHz as displayed in figures 
and ^ are performed in ~100 s and can be used only 
to characterise steady-state situations. A simple method 
to track the dynamics of buildup or decay of M in an OP 
experiment consists in actually repeating identical exper- 
iments with the probe laser frequency tuned to different 
transitions. The result of the analysis of two successively 
performed absorption recordings is given in figure The 
probed lines in these measurements are the two most in- 
tense lines in figure p8|a. The recorded signals are obtained 
with a standard modulation of the RF discharge excitation 
and a phase-sensitive analysis at /rf of the transmitted 
probe intensity. A higher nuclear polarisation is obtained 
here for two reasons: the gas pressure is lower, and a fibre 
amplifier is used to increase the pumping beam power . 

Small changes of the number density ns are derived 
from the recordings. They result in part from a slow drift 
in the RF discharge conditions, since in the analysis ^3(0) 
is obtained from the absorptions measured for the M=0 
situations at the beginning of the recordings. This drift 
accounts for the ^\ % difference observed over the total 



duration of the recordings (>700 s). Another part (~2 %) 
of the change of in figure ^ is clearly correlated to 
the value of the nuclear polarisation. It is known that the 
electronic orientation in the 2^S state has an influence on 
the rate of ionising collisions in the plasma: 

He* -I- He* ^ He + He+ + e" (55) 

(which would ultimately be forbidden in a fully polarised 
system |5^). Such an effect probably explains the system- 
atically observed variations of with M (here a small in- 
crease) . However these variations (even their sign) strongly 
depend on the RF discharge intensity |^ and are thus dif- 
ficult to quantitatively predict. This makes the simultane- 
ous measurements of M and 713 of special interest for an 
accurate characterisation of OP kinetics. 



3.3.3 Discussion 

The frequency scans in figures p8| a and p8[ b illustrate the 
difference arising from the chosen absorption measure- 
ment method: modulation of the RF discharge intensity 
and detection at /rf, or modulation of the probe fre- 
quency and detection at 2/probe- The signal-to-noise ratio 
of the phase-sensitive detection at /rf is observed to be 
higher when both amplitudes of the modulations are op- 
timally set. Detection at /rf may hence be preferred, but 
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this choice is not mandatory (a fair enough signal-to-noise 
ratio can be obtained, e.g. in figure |29| ). 

Whatever the chosen detection scheme, a modulation 
of the probe frequency may be applied to obtain a fre- 
quency error signal from a phase-sensitive detection at 
/probe- This allows to precisely adjust the probe frequency 
on the centre of isolated absorption lines, but is affected 
by the presence of any neighbouring lines (the effect de- 
pends on their distance and relative intensity). Whenever 
ill-resolved lines are probed to deduce relative populations 
of individual states, a careful analysis must be performed 
and appropriate corrections must be applied to the mea- 
sured signals. These corrections are usually smaller for 
signals detected at 2/probe when the distance to the near- 
est line is smaller than the Dopplcr width A. Conversely, 
smaller corrections are required for a detection at /r,f if 
this distance is larger than 1.3Z\, as is the case for the Cg 
doublet in Figure ^ (the exact crossover distance depends 
on the relative amplitude of the perturbing line) . This may 
seem surprising since ill-resolved lines have narrower fea- 
tures with a detection at 2/probc. However it is only the 
central part of the second derivative of a Gaussian profile 
which appears to be narrower than the profile itself. Side 
lobes actually extend beyond the Doppler width, and as 
a result there may be a significant crosstalk between ap- 
parently resolved lines. Indeed, the choice of a detection 
scheme impacts both on the data analysis and on the ac- 
curacy of the results only when measurements at two fixed 
frequencies are performed: for a complete frequency scan 
of a set of lines, integration of signals obtained by a de- 
tection at 2/probe can be used to exactly reconstruct the 
usual absorption spectra. 

We have tested several methods to obtain time-resolved 
results with dual-frequency operation of the probe laser. 
The one used for the recordings in figure |3^, which consists 
in repeating once every measurement, is very simple to im- 
plement but is time consuming. However, since the mea- 
sured values of the atomic density 713 and 714 usually vary 
very little with the polarisation, one may alternatively ne- 
glect this variation and on ly p erform one absorption mea- 
surement. Using equation |53| or [s^ with its left-hand site 
set to 1, the spin-temperature (and hence M) can then 
be obtained from this single absorption measurement. For 
the measurement of figure this approximation would 
have resulted in a maximum relative error of 2 % on M 
and 1 % on the decay rate Ti . 

For situations where very accurate measurements are 
required, or where the plasma conditions and hence 
and 714 may substantially change in the course of an ex- 
periment, alternating operation of the probe laser on two 
frequencies can be obtained by two techniques. A temper- 
ature change of the DBR diode can be used to induce the 
required frequency change (see section 3.1.1 ), with a typ- 
ical response time is of order 0.1 s. Following a frequency 
change of a few GHz, the frequency is actually sufficiently 
stable to provide a meaningful absorption signal after ~1 s 
of dead time. This means that the period of the probe fre- 
quency changes cannot be shorter that a few seconds, and 
that the technique is not suitable to measure fast pump- 



ing transients of the kind obtained using powerful pump- 
ing lasers . A change of the current in the DBR diode 
can also be used to induce the required frequency change. 
This is best done at a high enough frequency for the tem- 
perature of the device to remain constant (>1 kHz), and 
very fast OP transients can be monitored with this tech- 
nique. However, the current changes also drive significant 
power changes of the probe beam. The full elimination 
of this effect in the absorption measurements requires a 
careful signal acquisition and analysis, and the absolute 
accuracy of the results is somewhat reduced compared to 
the previous techniques. 

We finally come back to the initial assumption made 
in this section, namely that relaxation and OP processes 
do not affect the ME-ruled distribution of populations in 
the 2^S state, and discuss their perturbing influence on 
the optical absorption measurements. 

Relaxation of the spin variables in the 2'^S state[] would 
drive the spin distribution to its true thermal (Boltzmann) 
equilibrium in the applied field in the absence of ME col- 
lisions. Even for B=1.5 T and for electronic spins, this 
equilibrium corresponds to a small polarisation (<1 %) 
compared to M, so that the effect of relaxation is to re- 
duce the orientation in the probed 2'^S state. This leads 
to under-estimating the true value of the ground state nu- 
clear polarisation with absorption-based optical measure- 
ments. 

Conversely, OP with an intense pump laser results in 
an overpolarisation of the pumped 2"^S atoms with respect 
to the spin temperature population distribution (this is 
the driving process which actually transfers angular mo- 
mentum to the ground state reservoir of atoms through 
ME c ollisions). We have indeed demonstrated in section 
^.2.4| that even a moderate laser intensity can induce strong 



deviations from equilibrium of the population distribu- 
tions in high magnetic field. However this overpolarisation 
is a local effect which can simply be measured, and elim- 
inated by displacing the probe beam so as not to overlap 
the pump beam. Figure ^ clearly shows that blocking 
the pumping beam to let the nuclear polarisation decay 
does not change the measured relative populations nor the 
inferred value for M . 

Experimentally, the effects of relaxation and OP which 
we have observed at B=0.0925 T are similar to those ob- 
served in low field and do not significantly affect our 
optical measurement technique. This remains to be quan- 
titatively assessed in higher fields and at higher pressures, 
for which hyperfine decoupling and additional relaxation 
processes might decrease the efficacy of ME collisions to 
impose a single spin temperature for all populations. Still, 
based on our systematic study of OP effects up to 1.5 T, 
we believe that an absorption measurement technique can 
provide reliable results, especially if one probes the atomic 
states least sensitive to OP effects, and thus most effi- 



^ The main origin of this relaxation lies in the limited lifetime 
of the atoms in the 2^S state in the cell, due to diffusion to the 
wall (faster at low pressure) , ionising collisions (equation |55| 
more frequent at high densit y of He* atoms) , 3-body collisions 
(formation of He2 molecules |23|, more likely at high pressure). 
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ciently coupled to the ground state polarisation (the ^He 
2^S states A*^ and A'^, as shown in figure ||). 

4 Prospects 

From a theoretical point of view, we have developed and 
tested convenient tools to numerically obtain the eigen- 
states of the 2^S and 2^P levels in ^He and ^He and ex- 
plore the characteristic features of the 1083 nm transition 
in an arbitrary magnetic field. We have also extended the 
quantitative treatment of ME collisions by deriving rate 
equations for populations of all sublevels of the state 
which are valid for isotopic mixtures in an arbitrary mag- 
netic field. The use of an elementary OP model for the 
case M=0 led to quantitative predictions which have been 
accurately verified experimentally. A more elaborate OP 
model, similar to the one introduced for pure ^He in low 
field [[7|, can now be developed. This work is in progress, 
based on the rate equations describing the effects of ME 
collisions and absorption / emission of light which we have 
obtained and on an adequate phenomenological treatment 
of relaxation phenomena. 

Experimentally, we have demonstrated that the use of 
a DBR laser diode and of simple techniques based on stan- 
dard or Doppler-free absorption measurement can provide 
a wealth of spectroscopic information on the 1083 nm tran- 
sition in helium up to high magnetic fields. 

Saturated absorption experiments (or velocity-selective 
OP measurements) have been performed in this work only 
as a means to perform Doppler-free measurement of tran- 
sition frequencies, to check for the soundness of our exper- 
imental techniques. However the complex signals which we 
have obtained, combining narrow lines and broad pedestals, 
also provide quantitative information on coUisional redis- 
tribution of orientation between atoms of different velocity 
classes in the 2'^S state. These collisions are a key pro- 
cess when OP is performed at high pumping intensity and 
saturation effects are important jGQl- A detailed study is 
currently being made using independent pump and probe 
laser frequencies to fully characterise the coUisional trans- 
fer of orientation and the saturation effects. This work will 
also provide useful guidance for designing improved laser 
sources for high power OP. 

A detailed study of several systematic effects on ab- 
sorption measurements has allowed us to assess the accu- 
racy of relative population measurements in various sit- 
uations. It could be expected that hyperfine decoupling 
would hamper angular momentum transfer by ME colli- 
sions at very high field. We have demonstrated that ME 
collisions actually have a very different efficacy for cou- 
pling populations in the 2^S state to the nuclear polari- 
sation in the ground state depending on the value of the 
magnetic field, the helium isotope and the implied 2'^S 
state sublevel. Based on these techniques and results, and 
using the optical measurement methods described in the 
last section of this article, a systematic study of OP of 
helium in non-standard conditions (high magnetic field, 
high pressure) can now be performed. One of the moti- 
vations of this study is the observation, mentioned in the 



introduction, that higher nuclear polarisations may be ob- 
tained by applying a stronger field (0.1 T) than usual. We 
believe that a systematic study will help to elucidate the 
nature of relaxation processes in the plasma, and may lead 
to reduce their effect in situations of practical interest for 
applications requiring hyperpolarised '^He . 

We wish to thank the C.I.E.R.M. in the Kremlin-Bicetre hospi- 
tal for hosting the experiment. One of us (T.D) acknowledges 
support from the Ecole Normale Superieure. PICS and Polo- 
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Appendix: tables and matrices 

The tables and matrices used in the body of the text are 
collected in this Appendix. 

Fine-structure term parameters are listed in tables ^ 
and H: 





2^Po 


2^Pi 


2^P2 


2"P 


2"Po 


Eo 










2^Pi 





El 







2^P2 








E2 




2"Po 


Eo 











2^?! 





Ei 





Em 


2^P2 








E2 





2'P 





Em 





A 



Table 4. mj-independent matrix elements of the fine- 
structure term Hfs in the {\J;mj)} basis, used to compute 
the 2^P state structure. Eo, Ei, E2 are the true energies of 
the 2^Po,i,2 levels in zero B field. The last 4 lines of the table 
provide the expression of the effective Hamiltonian explicitly 
considering the 2^^? — 2^P mixing, with additional parameters 
A, Em and E[. 





''He (MHz) 


^He (MHz) 


reference 


E2 














El 


2291.175 


2291.926 






dl 


Eo 


29616.950 


29616.675 




51 




A 


6. 1423 X 10'' 


6.1431x10^ 




36 




Em 


-17035 


-17037 




36 




E'l 


2295.899 


2296.651 


see caption 



Table 5. Values of the energies for "^He and ^He are given 
with reference to that of the 2'^P2 level of each isotope. We use 
the most recent and accurate values whenever possible, and 
mass effects are evaluated as corrections [^,^. Ei and E'l 
are almost exactly 1 MHz larger, and Em 1000 times larger, 
than the values of reference |36|, probably due to misprints in 
the article. E[ is not a free parameter since the eigenvalue Ei 
must be obtained. It is given by E'l — Ei + Efj/ {A — Ei) . 
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The Zeeman Hamiltonian, including relativistic correc- 
tions of O (a'^) , consists of a linear term Hi and a quadratic 
term H2. The linear term contains a correction with re- 
spect to the simple expression in equations ^ and |^: 



Hi — Hz 



(Al) 



in which the tensor operator O*^^^ involves L [ po||3^ , ^ . 
Due to this additional term, Hi is more easily expressed 
in the coupled | J, mj) basis, as is the quadratic term H2- 
H1+H2 has non-zero elements only between sublevels of 
equal mj, which involve five different g-factors in the 2^^? 
state, and two in the state, all listed in table III of 
1^ . In the present work, this full expression is used only to 
evaluate the error resulting from neglected terms involving 
gx and H2. The values listed in table |^ which are used 
for g'g and g'j^ are taken from js^, limiting the number 
of significant digits in order to be consistent with having 
neglected g'^ . The mass-dependence of is also computed 
from results of ||3^, separating out the mass-dependent 
term: 



m 



2Fi { M, 

1 H — - - I 

^/6 V^^4 



+ 5gL, (A2) 



in which Mi stands for the mass of either isotope (z=3,4) 
and 5gL is the correction for *He given in the reference. 



2^S 
2^? 



9s 



9L (*He) 



9l 



^He 



2.002237319 
2.002238838 











0.999873626 0.999827935 



Table 6. Values of the g-factors used in the simplified form of 
the Zeeman Hamiltonian given in equations ^ and ^. The nu- 
clear term in equation ^ involves 57=2.31748x10"^. The value 
of the Bohr magneton in frequency units is /iB = 13996.242 
MHz/T H]. 





mL=l 


mL=0 


mL=-l 




1 








TT 





1 













1 



Table 7. Selection rules allow transitions only if angular mo- 
mentum is conserved, depending on the polarisation state of 
the light, e\ — a-i, a- or vr, and on m_L. 



The matrix Fq which represents the contact interaction 
term of the hyperfine Hamiltonian Hhfs (equation ^) in 
the decoupled {\ms : rnj)} basis and the values of the hy- 
perfine coupling constants Ag and for the 2'^S and 
states respectively are given in table pc 
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1:-|- 




-1:+ 


1:- 0:- -1 : - 


l:+ 


1/2 








Q:+ 








l/\/2 


-1:+ 






-1/2 


1/72 
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1/72 




-1/2 
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l/%/2 




-1:- 








1/2 



Table 8. Matrix elements of Fq/As (equation M). Values 
of the hyperfine coupling constants of equation pTare As — 
-4493.134 MHz ||^ and Ap=-4283.026 MHz, 0.019% lower 
than the value the corresponding parameter C of reference |^ 
to account for singlet admixture effects (see text). 



Computed line positions and intensities for B—Q which 
are represented in figure ^ are listed in table ^. The results 
for "^He slightly differ from previously published data 
due to the use of a more accurate Hamiltonian. 



The transition matrix elements T^'^- and Ty are ob- 
tained using the eigenstates of the 2*^8 and 2^^? states and 
the selection rule operator T{ex): 

Tl,tl, = \{ms\T{ex)\Z,)\' (A3) 
T,,^\{A,\T{ex)\B,)\\ (A4) 

with matrix elements of T{e\) given by: 

(ms \T{ex)\ m'g,mL) = 5{ms, m'g) x S{ex, m^) (A5) 

for ■^He, and by: 

{ms ■■ mi \T{ex)\ m'g, mL ■ m'l) = 

S{ms,m'g)6{mi,m'i) x Z:{ex,mL), (A6) 

for ■^He, with S{ex, mL) as given in table 0. 





e/h (GHz) 


Intensity 


Ci 





0.03566 


C2 


4.51159 


0.37482 


C3 


4.95920 


2 


C4 


5.18076 


1.29767 


Cs 


6.73970 


1.29767 


Ce 


11.25129 


0.29185 


C7 


11.92046 


0.03566 


Cg 


32.60452 


0.29185 


C9 


39.34422 


0.37482 


Da 


38.53362 


5/3 


Di 


40.82479 


1 


Do 


70.44174 


1/3 



Table 9. Optical transition energies and intensities for B=0. 
The energies e are referenced to that of the Ci transition, the 
intensities are the partial sums of the transition matrix ele- 
ments Tij for each line, and obey the total sum rule of equa- 
tion ^. 
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mL,ms 


J;mj 


1,1 


0,1 


-1,1 


1,0 


0,0 


-1,0 


1,-1 


0,-1 -1,-1 


2;2 


1 
















2;l 
2;0 

2;-2 




1/^2 


1/V6 


1/V2 


2/V6 


1/V2 


l/%/6 


1/^2 

1 


i;i 

1;0 

i;-i 

0;0 




-1/72 


-1/V2 
l/\/3 


1/V2 


-1/V3 


-l/^/2 


1/V2 
l/\/3 


1/V2 



Table 10. Matrix elements of the transformation operator P4 (in equation The inverse transformation is given by the 
transposed matrix : P^^ {i,j)=P4 {j, i) . 



The matrix elements (J, mj IP4I mi, mg) of the trans- 
formation operator P4 from the (L, S) to the (J) repre- 
sentation for the 2^P state of '''He are given in table . 
The additional term Hf^f"^ in equation || is given in ta- 
ble |ll|. It only involves 2 parameters d and e, similar to 
parameters D /2 and E/h defined by the Yale group |Q 
but with values adjusted to effectively take couplings to 
singlet states into account with this simple model. 



of mp- With different notations, these expressions would 
be identical to those given in reference pj] but for a mis- 
print (a missing factor 1/2 in the last line of their equation 
3). Here, both and are obtai ned f or po sitive sign 
in the second lines of equations AlC and All. The field 
dependence of all energies is shown in figure ^ in the main 
text. The corresponding eigenstates can be simply written 
introducing the field-dependent angles 9+ and 6^: 



Given the total Hamiltonian in the 2^S state, which 
couples only states \ms '■ mi) of given mp — ms + mi, a 
simple analytical expression of the eigenstates and of the 
energies can be derived. To simplify the following expres- 
sions, we introduce the reduced magnetic field b (note that 
< 0) and the nuclear correction e: 



b = -gsf^BB/As 

e = gif^g's- 



(A7) 
(A8) 



The mp = ±3/2 states, | — 1 : — ) and |1 : +) , are eigenval- 
ues with energies (1/2 ± /3 (1 + e)) Ag. In the mp = -1-1/2 
subspace, the matrix of the reduced Hamiltonian Hq/As 
on the basis {|0 : , |1 : — )} is: 



-be I/V2 
I/V2 -&(l-e)-l/2 



(A9) 



(the matrix for mp — — 1/2 on the basis { 1 : — ) , | — 1 ; -|-) } 



is obtained replacing /3 by -/3 in equation A9). The level 
energies are: 



h,l 



As 



4 

±79 



(1 + 26 



46 + 462-86e-1662e(l-e)) (AlO) 



h,l 



Ac 



(1 - 26 



± \/9 - 46 + 462 ^ _ i6^2£ (1 _ g))^ (All) 



in which the upper indices h and / refer to the higher and 
lower energies respectively, and the lower index is the sign 



tan 
tan 9^ 



A 



-V2{W'l/As + l/2 + b{l-e)) 

-V2{w'l/As~be) 

-sin 61+ |0 : +) +cose'+ |1 : -) 

cos6l+ |0 : +) +sin6l+ |1 : -) 

-sm9- |-1 : +) +cos6'_ |0 : -> 

008 6*- |-1 : +) +sin0_ |0 : -) • 



(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 



The variation of the angles (?_|_ and 9- in shown in fig- 
ure Bll The large amount of state mixing due to hyperfine 




Fig. 31. The level- mixing parameters sin 6+ and sin5_ are 
plotted as a function of the magnetic field B. Insert : asymp- 
totic 1/B decays at high field, with almost equal values of order 
1.1x10"^ at 10 Tesla. 
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mL,ms 


+ 






mL,ms:mi 


1 1 


0,1 


-1,1 


in n n 
i,U U,U 


1 n 
-1,0 


11 n 1 
1,-1 U,-l 


1 1 

-1,-1 




(i-\-ze 














U,i.-|- 






















—d+2e 


— 3e 








1,0:+ 




3e 




d 








0,0:+ 






-3e 







3e 




-1,0:+ 










-d 


-3e 




1,-1:+ 








3e 




d-2e 




0,-1:+ 










-3e 


4e 




-1,-1:+ 














-d-2e 


1,1:- 




V2(d+3e) 




-^/2e 








0,1:- 






V2(d-3e) 










-1,1:- 










-^/2e 






1,0:- 








V2d 








0,0:- 












2^e 




-1,0:- 














-\/2e 


1,-1:- 










6V2e 


V2(d-3e) 




0,-1:- 














V2(d+3e) 


-1,-1:- 

















Table 11. Matrix elements of Hj^Ys in equation |6[ Half of the matrix is explicitely given, the missing elements being obtained 
using the matrix symmetries with respect to the two diagonals (the same matrix element is obtained when all six signs of 
angular momentum projections are changed). The values of the parameters are d—— 14.507 MHz and e= 1.4861 MHz, slightly 
larger than the values of the corresponding parameters D/2 and E/5 ol reference (see text). 



coupling at low field (sin^ 9+ — 1/3, sin^ 6'_=2/3 for B=0) 
is strongly reduced above IT (sin^ 6'± ^ 0.012/5^). Our 
state naming convention in most of this article (Ai to Ag 
by increasing energies) introduces three field intervals due 
to two level crossings. The first crossing occurs at 0.1619 
T between |1 : +) and Afi, the second one at 4.716 T be- 
tween A^i and A*^. Altogether, three states are affected 
by these level crossings, as summarised in table 



B (T) 


< 0.1619 


0.1619-4.716 


> 4.716 


any B 


Al 








1-1:-) 


A2 








AL 


A3 




aV 


A'' 




A4 






AV 




As 


At 








Ae 




AX 



Table 12. Table of correspondence between the state names 
used in this Appendix and the energy- ordered naming conven- 
tion Al to As. 



ME matrices: we now provide the elements of the in- 
dependent ME matrices L, , , , E'^ and F'^ which 
appear in rate equations ^ to The six relative pop- 
ulations ai of the text are here designated by a , a5_, 

a[^_, a++, aL, a^-' This indexing convention is the same 
for the state names and the corresponding populations, so 
that table |l^ can be used to transform the indexing, or to 
determine the permutation of rows and columns required 
for B >0.1619 T between the matrices provided in tables 



^ to ^ and those appearing in equations |3^ to For 
more compact expressions, we use the following notations: 

c± = cos^ 6*^:, s± = sin^ 6'±. (A18) 

These matrices are not fully independent. In particular, 
the steady-state solutions of equations ^ and ^ given 
in equations ^ and are independent on the isotopic 
ratio R. Substitution of the steady state solutions yi = 
'^i^Gff.ttk of equation [33 in equation M can be used to 
show that E^ ^ E^ -W and = • G"*, two field- 
independent relations which can be directly verified on 
the matrices given in tables |l| to |l^. 



a — 


ai 




a++ 


ai 


4 


-1 


c_-s_ 


C+-S+ 


1 




S+-C+ 



Table 13. Matrix elements of operator L (equation 33|). 







ai 


1 

a+ 


a++ 


h 

a_ 


h 

a+ 




1 


C- 








S- 










S- 


c+ 





c_ 


s+ 


J/3 








s+ 


1 





c+ 



Table 14. Matrix elements of (equation ||). 
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2E^ 






yi 


ya 


2/3 




2/2 


y3 


a 


1 








-1 








a'_ 


c_ 


s_ 





c_ 


-s_ 





aX 





c+ 


s+ 





c+ 


-s+ 


a++ 








1 








1 


a'^ 


s_ 


c_ 





s_ 


-c_ 










s+ 


c+ 





s+ 


-c+ 



Table 15. Matrix elements of 2£"' and 2F'^ (twice those of 
and in equation BSl) . 





a — 




aV 


a++ 






a 


1 


C- 








s_ 





a' 


c_ 


„2 1 „2 
C_+S_ 


c+s_ 





2c_s_ 


s+s_ 


a\ 





C+S_ 


„2 1 „2 
C+ + S+ 


s+ 


c+c_ 


2c+s+ 


a++ 








s+ 


1 





c+ 


at 


S- 


2c_s_ 


C+C- 





ci+si 


c_s+ 







s+s_ 


2c+s+ 


c+ 


C-S+ 


4+4 



Table 16. Matrix elements of 2iJ^ (twice those of in equa- 
tion pSj). 





0-- 


a'_ 


4 


a++ 






a__ 


-1 


-c_ 








-s_ 







c_ 


„2 „2 
C_-S_ 


-C+S- 








-s+s_ 


«V 





C+S_ 


2 2 


-s+ 


c+c_ 





a++ 








s+ 


1 





c+ 




s_ 





-c+c_ 





2 2 


-c_s+ 












-c+ 




„2 „2 
S+-C+ 



Table 17. Matrix elements of 2_F^ (twice those of in equa- 
tion Iss] ). 
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